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Chapter  1  INTRODUCTION 


I.  Motivation 

Wave  propagation  in  random  media  has  been  studied  over 
the  last  two  decades.  Heightened  interest  in  this  problem  has 
been  mainly  due  to  the  large  number  of  problems  that  arise  in 
radio  physics*  acoustice,  plasma  studies  and  certain  other 
branches  of  physics.  From  the  physical  viewpoint  random  wave 
propagation  can  be  analyzed  at  two  levels t  macroscopic*  con¬ 
cerned  with  propagation  in  continuous  random  media  such  as 
turbulent  fluids  and  microscopic*  concerned  with  the  scatter¬ 
ing  of  waves  by  randomly  distributed  scatterors  such  as  elec¬ 
trons*  molecules*  rain*  blood  cells.  Wave  propagation  in 
continuous  random  media*  applies  to  such  problems  as  the 
scattering  of  sound  and  ultrasound  waves  in  sea  water*  light 
scattering  in  the  atmosphere*  radio  waves  scattering  in  the 
ionosphere  and  the  twinkling  of  stellar  images.  Wave  propa¬ 
gation  in  discrete  random  media  is  of  considerable  interest 
for  such  problems  as  molecular  scattering  of  light,  wave 
scattered  by  the  rain,  and  bioengneers  may  use  the  fluctua¬ 
tion  and  scattering  characteristics  of  a  sound  wave  as  a 
diagnoistic  tool.  The  abundance  and  variety  of  such  problems 
has  stimulated  development  and  refinement  of  statistical 
methods  for  calculating  wave  propagation  in  a  random  medium. 
This  paper  is  mainly  concerned  with  the  wave  propagation 
through  an  ionospheric  random  layer. 

When  a  random  wave  from  an  extraterrestrial  sources* 
such  as  a  ra<Jio  star  or  an  artificial  satellite*  passes 
through  the  ionosphere*  its  wavefront  and  amplitude  will  be 
distorted  by  the  density  fluctuation  in  the  ionosphere  plas¬ 
ma.  On  propagating  to  the  ground,  the  resulting  phase  varia¬ 
tions  cause  interference  to  occur  and  a  diffraction  pattern 
is  set  up  across  the  ground.  The  resulting  random  variations 


of  the  electromagnetic  field  are  termed  scintillations.  In 
consequence ,  satellite  based  communications  eudiibit  serious 
performance  degradation,  the  degradation  being  most  serious 
for  propagation  paths  which  transit  the  auroral  and  equatorial 
ionosphere.  The  importance  of  a  through  investigation  of  this 
problem  is  obvious. 

II.  General  Introduction  and  Review 

Propagation  models  are  required  to  provide  satellite 
communication  system  and  radar  system  designers  with  a  means 
to  translate  the  available  ionospheric  statistics  into  cons¬ 
traints  for  the  specific  systems  that  are  being  designed. 
Approximate  models  have  been  proposed  that  apply  under  res¬ 
trictive  conditions.  Analyses  of  propagation  effects  general¬ 
ly  start  with  the  scalar  wave  equation.  If  the  dominant 
scattering  irregularities  have  dimensions  much  larger  than  a 
wavelength,  depolaization  effect  can  be  neglected  (1,  2). 

Let  the  wave  field  be  represented  by  E  =  E  exp  (-iet),  where 
m  a  2/1  f,  f  is  the  carrier  frequency,  t  is  the  time,  E  is  the 
electric  field  vector  and  u  is  a  component  of  E.  If  the 
dielectric  properties  of  the  medium  change  slowly  in  time  in 
comparison  with  1/f  and  slowly  in  space  in  comparison  with 
the  wavelength  A  ,  then  with  a  monochromatic  source,  the 
propagation  of  the  wave  through  a  random  irregularity  region 
is  governed  by  the  scalar  wave  equation 

A  ^  a  A  A 

r*u  +  ft  £  u  -  O 

(1.1) 

A  A 

where  k  »  sS.  ,  £  is  the  dielectric  constant  and  both  U  and 

A  A  A 

£  are  random  variables.  £  enters  as  a  coefficient  of  the 
unknown  wave  function  U  •  This  is  the  root  of  all  the  mat he- 


matical  difficulties  of  the  theory*  since  we  don't  know  how  to 
find  an  exact  solution  of  such  a  wave  equation*  It  is  necessa¬ 
ry  to  apply  certain  approximations  that  make  use  of  the  small 
parameters  such  as  the  fluctuation  of  £  *  i.e.  the  deviation 
from  the  mean  value.  It  can  also  be  the  smallness  of  the  wave¬ 
length  in  comparison  with  the  dimension  of  the  inhomogeneities. 
In  practice*  two  types  of  problems  arise;  the  direct  problem* 
in  which  one  has  to  find  the  statistics  of  waves  propagating 
in  the  medium  from  the  known  statistics  of  the  medium,  and  the 
inverse  problem  wherein  one  draws  conclusions  about  the  proper¬ 
ties  of  random  inhomogeneities  from  the  measured  moments  of  the 
field.  Actually*  these  two  problems  cure  equivalent,  .i.e.  one 
needs  to  find  the  relation  between  the  statistice  of  the  medium 
and  the  wave  field. 

We  shall  assume  for  simplicity  that  the  medium  is  on  the 
average  homogeneous  and  stationary.  Removing  these  two  res¬ 
trictions  does  not  give  rise  to  any  difficulties  in  principle. 
Under  conditions  usually  encountered  in  ionospheric  propaga¬ 
tion*  the  energy  is  scattered  into  the  forward  direction  and 
(1.1)  can  be  replaced  by  the  parabolic  equation! 

+  Vf*U  +  i  0  =  O 

(1.2) 

where  k»  «  ka<€>  *  u  ■  U  exp(ik*z),  z  is  the  direction 

of  propagation*  and  Vf1  ■  +  *  the  transverse 

Laplacian.  Both  the  wave  equation  (1.1)  and  the  parabolic 
equation  (1.2)  are  stocdiastic.  The  wave  properties  of 
interest  are  the  average*  variance*  and  higher  moment  of  U  . 

If  B  and  U  were  uncorrelated*  these  equations  could  be 


solved  easily.  However*  £  and  U  are  generally  correlated, 
the  approximations  used  to  provide  solutions  to  (1.1)  and  (1.2) 
are  all  directed  towards  modeling  the  statistical  relationship 

rnm  A  mm 

between  £  and  U  .  If  £  is  small  enough,  then  we  can 
naturally  resort  to  the  method  of  perturbations,  and  expand 
in  a  power  series  in  £  ,  or  more  exactly  in  (€*)^  . 

A  pertubation  series  solution  to  (1.2)  can  be  formed  that 
separates  £  and  U  where  U  describes  the  fluctuation  in 

A  A  A  mm 

tJ  ,  (7  »  <U>^+  V  .  The  first  term  in  the  perturbation  series 

for  V  is  the  Born  or  single -scattering  approximation,  which 
applies  only  when  both  £  and  U  are  small.  When  U  is  not 
small,  a  large  number  of  terms  in  the  perturbation  series  must 
be  summed.  The  nth  term  of  the  series  describes  the  n-fold 

^  mm 

scattering  and  contains  the  n-fold  product  £  (ri ) . . .  £ (rn ) . 

,  A 

Thus,  in  calculating  the  average  <U>,  we  have  to  know  the 

mm  mm  mm 

moments  <  £  (rt)...  £  (r*)>  of  £  of  all  orders.  The  reno- 
malization  schemes  of  multiple-scattering  theory  attempt  to 
solve  this  problem  by  a  selective  summation  technique  that 
leads  to  readily  evaluated  expressions  for  the  moment  of  U 
of  interest.  v 

The  Born  approximation  is  shown  to  be  valid  when  1c,  1  is 
small,  where  1  is  the  characteristic  scale  of  the  turbulent 
medium.  As  k0  1  increases,  we  must  either  take  into  account 
the  higher  order  terms  in  the  perturbation  series  or  go  over 
to  other  approximate  methods  which  deal  with  multiple  scatter- 

A 

ing  to  some  extent.  Rytov  proposed  that  an  equation  for  4  ** 

A  A 

in  U  be  used  in  place  of  (1.2).  The  equation  for  W  is 
then  given  by 
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A  00 

the  equation  has  seperated  *±'  and  £  so  that  the  correlation 

A  M 

between  end  £  is  not  explicit*  A  perturbation  series 
expansion  for  ,  representing  the  fluctuation  in  'i'  ,  can 
be  formed  to  solve  the  nonlinear  equation  (  in  ).  The 
zeroth  order  equation  represents  free  space  propagation.  The 
1st  order  equation  for  $  can  be  solved  yielding  the  Rytov 
approximation ,  which  is  also  called  the  method  of  smooth  per¬ 
turbations.  The  Rytov  approximation  is  a  weak  scintillation 
approximation  that  hold  for  small  <X*>  ,  where  X  is  the  flu¬ 
ctuation  of  the  logarithm  of  the  amplitude,  it  has  a  wider 

range  of  validity  than  the  Born  approximation. 

« 

The  generally  accepted  model  for  ionospheric  propagation 
thru  a  turbulent  ionosphere  has  been  the  thin  phase  diffrac¬ 
tion  screen  model  (3).  The  ionospheric  irregularities  perturb 
the  phase  of  the  field  at  the  layer  and  their  effect  upon  the 
wave  beyond  the  layer  can  be  computed  using  diffraction  theory. 
Mercier  (4)  and  Briggs  and  Parkin  (5)  introduced  a  Gaussian 
correlation  function  to  describe  the  phase  fluctuation  at  the 
screen  and  related  all  higher  moments  to  the  second  moment. 
Spectral  analyses  of  observed  phase  and  amplitude  fluctuations 
have  shown  that  the  electron-density  irregularities  have  a 
power  spectrum  that  may  be  characterized  by  a  power  law  shape. 
Rufenach  (6)  extended  tha  phase-screen  theory  using  a  power 
law  spectral  shape  for  the  irregularities. 

For  the  weak  scintillation  case,  the  Born,  Rytov  and  the 
thin  phase  screen  approximations  are  all  applicable.  For 
strong  scintillation  and  a  thick  layer,  the  effect  of  multiple 
scattering  on  transionospheric  signals  has  to  be  taken  into 
account. 

Several  techniques  have  been  proposed  to  deal  with  the 
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strong  fluctuation  problem.  They  are  the  diagram  method 
(DeWolf  1968,  1977;  Tatarski  1971;  Frisch  1968)  (1,  7,  8)  the 
integral  equation  method  (Brovn  1971,  1972),  (9,10),  including 
the  Dyson  and  the  Bethe-Salpeter  equations,  the  extended 
Huygens-Fresnel  principle  (Kon  1970;  Clifford  1974)  (11,  12), 
and  the  parabolic  equation  method  (Tatarski  1971;  Furustsu 
1972)  (1,  13).  In  this  paper,  we  shall  apply  the  parabolic 
equation  method  to  the  case  of  radio  wave  on  propagation 
through  an  ionospheric  slab.  First  the  statistical  moments  of 
waves  are  obtained  inside  the  random  slab.  Then  the  free  space 
propagation  of  resulting  randomly  modulated  waves  is  analyzed 
from  the  bottom  of  the  slab  to  a  ground  receiver.  The  geometry 
of  this  problem  is  shown  in  Fig.  1. 


Chapter  2  THEORETICAL  ANALYSIS 


I.  Problem  Description 

Ionospheric  scintillation  has  been  observed  at  frequen¬ 
cies  ranging  from  10  MHz  to  6  GHz.  The  irregularities  in  the 
ionosphere  which  produce  the  fluctuations  are  believed  to  be 
mainly  in  the  F-regiont  the  layer  of  strong  irregularities  is 
often  100-500  km  thick  ranging  from  200-700  km  in  altitude* 
The  electron  density  fluctuations  are  often  of  the  order  of  a 
few  tens  of  per  cent  and  can  be  as  high  as  70  per  cent.  The 
F-region  irregularities  are  usually  regarded  as  a  stochastic 
process  and  characterized  by  three  dimensional  random  func¬ 
tions  which  are  assumed  spatially  and  temporally  stationary 
over  the  increments  of  interest.  One  defines  the  spatial 
correlation  function  BM(r")  of  the  density  function  by 

<N<n> 

N  ~  <  N  1  (2.1) 


where  r  *  (x,y,z)  is  the  position  vector,  r**  ■  r.  -  r.  is  the 
correlation  lag  vector,  <N  ”is  the  rms  density  fluctuation, 
and  <  >  denotes  the  ensemble  average.  The  relationship 
between  the  three  dimensional  power  spectrum  of  the  density 
fluctuation  ^(k)  and  the  density  correlation  function  is 
given  by  the  Fourier  transform  palrt 
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^  ^  =  Hi 


(2*3 ) 


the  pover  spectrum  is  often  characterized  by  a  vide  range  of 
wave  numbers,  and  hence  designated  as  a  vide -band  process, 
whereas  a  narrow  range  of  wave  numbers  associated  with  a  Gau¬ 
ssian  spectrum  is  designated  narrow-band  process.  The  wide¬ 
band  process  is  important  since  it  is  inferred  from  the  in- 
si  tu  results.  In  the  F-region,  the  vide  range  of  wavenumbers 
corresponds  to  dimensions  ranging  from  a  fev  meters  to  hun¬ 
dreds  of  kilometers .  if  a  Gaussian  spectral  shape  is  assumed, 
it  is  characterized  by  one  dominant  scale  approximately  equal 
to  the  Fresuel  wave  number,  (about  2km”1)  (Brigg  and  Parkin 
1963). 

Refer ing  to  Fig.  1,  let  us  consider  a  time  harmonic 
radio  wave  incident  on  a  region  of  ionospheric  irregularities 
at  z  *  z0  .  L  is  the  thickness  of  the  slab.  Inside  the 
irregular  region,  the  relative  dielectric  permittivity  is 
given  by 


i(Y)  =  £.<?)[!  +  S(  r>J 


(2.4) 


where  one  models  a  collisionless  ionospheric  plasma  by 


£. (})  s  1  - 


ui.l  a) 


(2.5) 
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1  )  is  the  plasma  frequency  for  the  background  electron 
density  profile,  N(r)/N0  is  the  percentage  fluctuations  of 
the  electron  density.  In  the  following,  we  shall  assume  that 
N(r)  is  a  homogeneous ,  isotropic  random  field  while  is 
taken  as  constant. 

The  wave  first  travels  through  the  ionospheric  slab  of 
length  L,  at  z  »  z^  the  modulated  wave  enters  the  free  region, 
and  it  is  detected  at  z  =  z  ^  •  Therefore,  there  are  two 
regionsi  in  the  first  region,  the  incident  wave  at  z9 is  deter¬ 
ministic  and  the  outgoing  wave  at  zt  becomes  stochastic.  In 
the  second  region,  the  entering  wave  at  z4  has  a  stochastic 
modulation  but  the  medium  is  deterministicv 

Writing  the  scalar  wave  field  as 


E 


U  <r) 

00 


e 


(2.7) 


where  the  exp(-iwt)  dependence  is  omitted,  k*6# —*  k#  ,  and 
u(r)  is  the  complex  amplitude  of  the  wave.  Substituting  (2.7) 
into  the  wave  equation  (1.1),  we  obtain  an  equation  for  u  (r) 


-  h  U  tn 
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+  V*U  +  *:Z(l>U  =  O 


(2.8) 


where  X7p  z.  V  —  g“jT  .  If  the  complex  amplitude  u  varies 
markedly  over  distances  of  the  order  of  the  inhomogeneity  scale 
1,  the  second  derivative  3  U  /  9|l  is  of  the  order  of  u/12. 
On  the  other  hand,  the  term  2ilc03u/a|1  in  (2.8)  is  of  the 
order  of  u/Al.  Therefore,  for  X  <<  1,  the  term  9  W-/32*  is 
small  compared  to  the  first  term  in  (2.8).  Thus,  one  can  re¬ 
place  Vau  in  (2.8)  by  the  transverse  Laplacian  u  and 

obtain  the  parabolic  equation 


+  VflU  +  6.zZ(r)U  = 


(2.9) 


This  is  the  starting  point  of  our  analysis,  in  the  later 
sections,  we  shall  examine  the  statistical  properties  of  u  j 
they  are  specified  by  the  infinite  set  of  correlation  functions 


=  <u  -  w  (.i.U)  w*(3»r.V  -u*  »•*•')> (2-10> 


The  general  theory  of  multiple  scattering  has  been  deve¬ 
loped  by  many  authors.  (Bourret  1961,  Apresyan  1973,  Tatarski 
1971,  Marcuvitz  1974,  DeWolf  1979).  The  Greenes  function  me¬ 
thod,  which  had  previously  been  investigated  in  quantum  field 
theory,  was  applied  to  the  problem  of  wave  propagation  in  a 
random  medium.  People  have  used  the  Dyson  equation  for  the 
average  field,  and  the  Bethe-Salpeter  equation  for  the  cova¬ 
riance  B  *(u(r^ )  u(r»)^.  However,  one  cannot  derive  such 
closed  equations  by  averaging  the  original  differential  equa¬ 
tions  for  the  field  because  moments  of  different  orders  are 
coupled  together.  Bourret  (14)  was  the  first  person  to  apply 
the  graph  technique  to  scattering  of  waves  in  a  continuous 
fluctuating  medium.  He  assumed  that  the  parameters  of  the 
medium  fluctuate  according  to  a  normal  law  and  are  statistical¬ 
ly  independent  of  the  sought  field.  This  led  to  an  approximate 
expression  for  the  effective  average  dielectric  constant  pro- 
portional  to  the  correlation  function  <  £  (r*  )  £  (r» )  y  of  the 
medium.  A  more  detailed  derivation  of  the  equations  for  the 
average  and  2  point  correlation  assuming  a  normal  law  for  a 
fluctuating  medium  is  due  to  Tatarski,  and  to  Frisch  who 
admitted  deviations  from  the  normal  law. 

The  utility  of  series  representations  is  dependent  on 
the  rapidity  of  convergence.  Operator  techniques  employed  in 
the  spirit  of  the  formal  theory  of  scattering  was  shown  by 
Marcuvitz  (1974),  to  provide  alternative  and  formally  exact 
representations  of  turbulent  field  quantities  and  their  n- 
point  ensemble  averages.  In  appropriate  parametric  ranges 
these  formally  exact  expressions  are  expandable  into  rapidly 
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convergent  series.  In  this  section,  we  shall  apply  this 
technique  to  the  average  equation  and  investigate  its  conver 
gence. 


In  abstract  the  wave  equation  (2.9)  can  be  written  as 


(L.-V)f  =  0 


(2.11) 


re,  in  a  f  ,  z  space  the  unperturbed  operator  L«  .  and 
perturbation  V  ,  whose  average  4.  V  >  *  0  ,  are  represented 


(2.12) 


equation  (7)  is  a  generic  operator  equation  for  the  field,  vith 
the  initial  condition  acting  as  a  impulsive  source  applied  at 
z  *  0.  It  is  convenient  to  define  a  stochastic  Green's  func¬ 
tion  such  that 

(  Lo  ~  W  G*  ^  1  (2.13) 

where  in  f  ,  z  space 

Qf  - ►  J  f  ( t,i;  J3 

'  i  - *•  t-n  b  a-}')  <*•“> 

0+ 
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and  &  (x)  is  the  Dirac  delta  function.  For  uniqueness,  6  will 
be  distinguished  by  the  requirement 


§  (f,>,  f  ,  V)  =  o 


(2.15) 


If  we  introduce  G0  as  unperturbed  operator  inverse  to  the 
operator  L0  ,  and  defined  by 


L  o  Cr 


0 


1 


(2.16) 


where  the  domain  of  L0  is  such  that 


a. - *  s  o 

•  V 


(2.17) 


A 

as  is  the  case  for  6. 

A 

The  ensemble  average  G  of  G  is  taken  to  satisfy  v 

(  L.O  “  Vc  )  Cf  =  1  (2.18) 

and  depends  on  the  "smoothed"  scattering  operator  Vc  defined 
below.  Eq.  (2.18)  provides  a  nonlinear  defining  equation  for 
6,  the  nonlinearity  arising  from  the  dependence  of  Vc  on  G. 


ytiaMMii 
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A  stochastic  operator  Tc  ,  with  zero  ensemble  average, 
is  defined  by 


G  :  G  +  G  Tc  G  (2-19) 


which  represents  the  multiple  scattering  in  the  smoothed  back- 
ground  G.  One  can  derive  the  expression  for  Tc  in  terms  of 
G  and  a  power  series  in  V  to  4th  order  in  V  t  (Appendix  2) 


s  y  +  v  Gr  V  ±  V  Gr  V  Gf  V 

-  <vGrVQV>  -<V<=tV>QV 
-vGr<vQv)>  +  v  G  V  V  Gr  V 

7Qv>-vQ<vQv>Gv 
t  <vGr<v$V>GrV  >-<vCrV>GrVGrV  (2*20) 

-  vGv  Q  <v(kv  >  t  2  <  vfiv?  G<vGv> 

- <vGvqv>Gv -vGKvGvGv  > 

+ . 


and  correspondingly  the  operator  Vc  =  <  VGTC  /  is  represented 
as 


vc  =  <vQV'>t<VrGV<TV’>t<V'(fl/(TV&V> 

-  <\7G<?<kv>(iv> 

-  <  VG  v>  §  <  vG  v>  + 


(2.21) 
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the  expansions  are  in  terms  of  the  average  operator  G,  which 
in  contrast  to  G 0  is  nonsingular.  They  are  more  rapidly  con¬ 
vergent  than  the  perturbative  expansions  in  terms  of  G0. 

As  a  simple  example,  used  to  study  wave  propagation  in 
a  one  dimensional  random  medium,  let  us  consider  a  differen¬ 
tial  equation  of  the  following  typet 


a  d 


(2.22) 


If  V(z)  =  V-£  (1  +W(z)),  where  is  a  real  mean  frequency 
and  P  (z)  is  a  centered  stationary  and  Gaussian  random  func¬ 
tion  of  z,  then  it  represents  a  randomly  modulated  oscillator 
whose  correlation  function  is 


B(3-3’)  =  rm> 


(2.23) 


The  model  will  be  used  as  a  check  for  the  approximate  expan¬ 
sion  in  equation  (2.21). 

The  Generic  stochastic  Green's  function  for  (2.22)  is 
defined  by 


d  }  (2.24) 
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To  find  the  average  Green's  function,  we  must  calculate 


Gtd.}') 


where  U(u)  is  the  step  function  to  insure  the  uniqueness  of  G. 

Here  C^V(Z)dT  ,  being  a  linear  functional  of  the  centered 

«• * 

Gaussian  random  function  V(z),  is  a  centered  random  variable 
y  ,  hence 


f<yJ> 


(2.26) 


Since 


<  (  £  v(«Jr)*> 

P  P  BU-Odrdr' 

J3'  -'3' 


(2.27) 


We  obtain  the  final  result 


Or  c  5.3')  =  -i  e 


l 


U  (  )  (2.28) 
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If  we  use  the  correlation  function 


=  A  b  (3-3') 


(2.29) 


the  integration  in  Eq.  (2.28)  can  be  performed  exactly  to  give 


-4<a-r; 

G(j.j')*-ie.  ua-}') 


(2.30) 


If  we  approximate  VCG  by  taking  the  1st  term  in  eq.  (2.21), 
i.e. 


VcG  =  <v&v>  <r 


(2.31) 


or 


Vc(r  = 


•J 


(2.32) 


Substitute  eq.  (2.30)  into  (2.32),  one  obtains 


4  GrU.i)  GU.V)V U-i') 

2  (2.33) 


VcG  = 


18 


since 

Cr  (  3/  5)  =  “ 


(2.34) 


therefore 


VcG  =-fAG<}.3'jUC3-5'j 


(2.35) 


the  solution  of  eq.  (2.24)  can  be  easily  obtained  to  give 

-4  <3-3') 

G-u.n=-;e  u(3-n  ,2 


which  is  identical  with  the  exact  solution  (2.30).  It  is  an 
interesting  result  and  leads  to  the  conclusion  that  t  If  V 
is  a  centered  stationary  and  Gaussian  random  function*  with 
a  correlation  function  defined  by  (2.29)*  the  first  term  in 
eq.  (2.21)  gives  the  exact  expression  for  V^G  .  i.e.  the 
higher  moments  terms  do  not  contribute  to  the  expansion. 
Furthermore*  since 


<vG'V>z<v<*.V>  +  <vGfoV(k0VG.V> 

+ . 


(2.37) 
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it  can  be  shown  that  in  this  particular  case 


<V«rV>  *.<VS.V> 


(2.38) 


Let  us  consider  a  mo  re  complicate  correlation  function! 


-iLli 

B(3-3')=  T 


(2.39) 


the  integration  in  eq.  (2.28)  can  be  evaluated  exactly  to  give 


(2.40) 


for  short-range  correlation. 


a  -i'T  <3-3')  (2-41’ 

V* 
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Therefore 


-Sril-t) 


(2.42) 


One  can  examine  the  1st  approximation  in  eq.  (2.31)  by  substi¬ 
tuting  (2.42)  into  the  expression  for  VCG: 

VcQ-  (3)  vd)>  G  fS.l;  Gf 

U(3-i)U  (3-3'Jdf  g-?- 

=  -  6iTe'<lTt,'s',f  l- e  ”r  J  (2'<,3) 

while  the  exact  solution  for  G  gives  rise  to  the  expression 

VCG  =  ijj  -  * 


=  - e 


(2.42) 


Compare  (2.43)  to  (2.44),  one  observes  that  VCG  *  ( VGV  )is  a 
good  approximation  when  the  correlation  range  is  small. 

Thus,  one  concludes  that  for  a  centered  Gaussian  random 
variable  V  with  a  delta  correlated  2-point  correlation  func¬ 
tion,  VcG  can  be  exactly  expressed  by  the  1st  term  of  the 
series  expansion  (2.21),  i.e.  VCG  *  <  vgv>  G  »  (  VG#V)G.  For 


■£pr»l4ptK.l.'  9t 


S'K 


nail  range  correlation  function,  ».<VGV>  G  is  a  good 
approximation.  The  representation  (2.31)  can  also  be  applied 
to  the  case  of  large  scale  correlation,  but  then  the  fluctua¬ 
tion  oust  be  small,  otherwise ,  the  higher  moment  terms 
have  to  be  summed  to  get  a  good  expansion. 

111.  Markov  Approximation 

In  order  to  proceed  further,  additional  approximations 
for  the  random  medium  statistics  must  be  considered.  A  very- 
useful  approximation  is  the  Markov  approximation  which  gives 
correct  results  in  cases  of  interest  even  in  the  limit  of 
strong  scintillations.  The  mathematical  representation  and 
physical  interpretation  of  the  Markov  approximation  is  based 
on  the  following  observation! the  correlation  of  the  dielectric 
constant  in  the  transverse  direction  £  has  a  direct  bearing 
on  the  transverse  correlation  of  the  field  but  the  correlation 
of  the  dielectric  constant  in  the  direction  of  the  wave  pro¬ 
pagation  has  little  effect  on  the  fluctuation  characteristic 
of  the  wave.  Hence,  one  can  assume  that  £  (r)  is  delta  corre¬ 
lated  in  the  direction  of  propagation,  v 

<e  (j.f)  £(?'£')>=  &  <}-}')  A  (l-t') 


i.e.  one  can  treat  wave  propagation  in  a  random  medium  as  a 
Markov  random  process.  It  is  equivalent  to  assume  that  the 
turbulent  eddies  are  like  flat  disks  oriented  normally  to  the 
propagation  path.  Tatarski  (1971)  gives  a  derivation  of  this 
approximation  which  provides  the  limits  of  applicability  of 
the  Markov  approximation.  For  typical  ionospheric  parameters. 
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these  conditions  are  generally  satisfied  for  an  incident  wave 
with  frequency  about  twice  the  maximum  plasma  frequency  or 
higher. 

The  relationship  between  the  function  A(  9  -  f  )  and  the 
power  spectrum  $  (k)  is  given  by  the  following! 

f  ,  ,  U-(P-f') 

A(f-f’)=  in  \$(i)e  '  n 

~  ~  J  *  (2.46) 


If  the  turbulence  is  isotropic,  i.e.  $  (k)  °  <§  (k)f  then 
get 


we 


A  (P-f')  =  (2nf  j  J0 

'■'O  (2. 


47) 


A( P  )  is  also  related  to  the  correlation  vf unction  B(  z, f  ) 
through 

A(r-f)=/'B<  »■£-£')« 

The  limit  of  application  of  the  Markov  approximation  was  inves¬ 
tigated  by  Tatarski  (1971)  via  variational  techniques.  We 
derived  identical  results  using  the  operator  method  as  present 
in  Appendix  3. 


. 
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IV.  Equation  for  the  Mean  Field 

In  summary,  our  analysis  of  strong  scintillations  will 
based  on  the  following  assumptions: 

(1)  The  parabolic  approximation,  (limits  of  applicability  of 
the  parabolic  approximations  will  be  derive  in  Appendix  4). 

(2 }  £  (r)  is  delta  correlated  in  the  direction  of  propagation. 

Let  us  take  the  ensemble  average  of  parabolic  equation 

(2.9) 

vf2<<?>  +  *.2<£' P>=  o 

(2.49) 

will  be  expressed  as  follows: 

=  vc(r)<^r)> 


. . „  3<^>  , 

2t.fi.-r-: —  + 
o  i 


The  last  term  in  (2.49) 


<£(  r)  *p<r)> 


which  is  valid  if  the  incident  field  at  z  *  0  is  deterministic. 
One  notes  that  £  (r)  ■  2n(r)  is  the  fluctuation  of  the  refrac- 
tive  index.  To  get  a  differential  equation  for  ^  ^  (r )^,  under 
the  assumptions  (1)  -  (2),  one  obtains  the  following  expre¬ 
ssion  for  Vci 


Vc  *  <  vGt  V>  =  <  vQeV> 


(2.50) 


.  ^  H  •  raMawlillH 


If  the  fluctuation  £  (r)  is  a  centered  Gaussian  random  field 

00 

and  characterized  completely  by  the  correlation  function 
<£(ri)£(rj)^  #  then  (2.50)  is  exact. 

To  ascertain  Vc,  the  unperturbed  operator  G0must  be  found. 
Therefore,  ve  have  to  find  the  Green's  function  for  the  homo- 
geneous  equation 


•  A  +  j.  j!i_  -  o 
A  +  2*. 


(2.51) 


Setting 


a  rT~ »  4  , , 

t  (f,3)  =  f  a,i)  e  - '  J i 

—  oo 


(2.52) 


One  can  write  equation  (2.51)  as 


•  iilM)  _  V  £  (tj)  =  o 


si 


(2.53) 


From  (2.51),  ve  define  G#  by 


•  3  0-.  ,  1  9*  Go 

*7T  +  2*.  3P* 


=  £  (  r-r’;  u.i 

^  V  v 


Setting 


G0ru,f  :j> 


iA-(f-f') 

rj_& 


f2/t)x 


where 


Ga(f-5J',3’)  »  o 


3  <  3' 


From  (2.55),  we  obtain 


xili  -i! 

a  }  2  fe< 


9*  =  S  Ci-3') 


=  0 


3  <»3 


For  z>z*,  we  note  that 


%  IIS')  = 


Using  (2.57)  and  (2.58),  one  sees  that 

%(}.}')  -  -  *  as  5-?' 


(2.55) 

(2.56) 

(2.57) 

(2.58) 

(2.59) 

(2.60) 


the  above  is  also  a  good  approximation  provided  k2  |z-z'|  /2ke«l . 
Therefore,  by  using  (2.50),  we  obtain 

-  -i$02  f*  <n 

°  A-oo  ~  (2.63) 

where  n  *  £/2.  Setting  z  -  z*  «  ^  ,  we  obtain 

Vc*-it>.2  fa<n(f.j;7l(f,3-t)>J7 

*  °  (2.64) 


Applying  the  Markov  approximation,  one  obtains 
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vc  =  -i 4.2  f~  Sen  A  (3.0)  c/7 


Therefore,  the  average  equation  (2.49)  becomes 

;  3Cf>  .  _L_  +  A  <b°XH>>=  O 

+  2&.  9f*  2  , 


I 


\ 


V.  Parabolic  Equation  for  the  Mutual  Coherence  Function 


The  parabolic  equation  (2.9)  and  its  conjugate  can  be 
written  asi 


(2.67) 


+  (2-68’ 


Multiplying  (2.67)  byt|/(z,f')  and  (2.68)  by  )  and  sub- 

*0  ** 

tracting  the  second  equation  from  the  first,  we  obtain* 
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*  A  * 

5=  o 


I 
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The  mutual  coherence  function, <  $  (z,  f  )  (z,  f '  )  >  is 

^  v  ^ 

defined  as  the  cross-correlation  function  of  the  fields  V  and 

A  . 

if>  m  a  direction  transverse  to  the  direction  of  propagation. 

Using  the^  operator  notation  described  in  section  II,  we 
rewrite  M  =  f  •?'*  in  the  forms 


(  L.  -  V  )  M  =  0 


(2.70) 


where 

^n  +  A(vr2-vt‘) 

<}./)'  n(}. n] 

V 

An  unperturbed  Green's  function  may  be  defined  as 

+  jkJvrvr)] 

=  i>  Ct-f)  i  (f'-F')S  Ci-%) 

~  •*  **  ** 


(2.72) 


Applying  the  same  analysis  as  described  in  section  IV,  one 
finds  approximately  that 


QAf.f',  lili',1)  = 


SfM)i(f-f') 


(2.73) 


The  smoothed  scattering  operator  for  the  ensemble  averaged 
Green's  function,  defined  as  in  (2.38)  can  be  expressed  by 


Vc  =  <  v  QV>  -  <v$.V> 


(2.74) 


and  written  as* 


Vc  s  - i *.*  jTf %  <  C  "  < »  - "  < l'  Wj  H l -f) 

-oo  "'00 

i  ( t-  VI n  (f ; S) ■ - ■ n  (l \  SjJ >  dfjf 'd  * 


Setting  z  -  z  *  *[  ,  one  obtains 

vc » -iVjyscvAM-swAfhf-rWl  w 
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Applying  the  Markov  approximation,  one  obtains 


[Albo)-A(i.r-i')] 

L  U  (2.77) 


Equation  (2.70)  then  becomes 


{*A  +4(V'  V) + 

s  O 


where 


A  ~  •  (2. 
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VI.  The  Equation  for  Fourth  Order  Moment 

A 

On  considering  the  parabolic  equation  for  if*  (z,  ?,  )  , 

A  A  ^ 

<z,fs),  (J/  (z,f4),  respectively,  we  can  obtain  an 
^  « 
equation  for 

m  =  ifa.w'P 

(2.80) 


as  in  the  form 


[  u  -  V  ]  ^4  -  O 


(2.81) 


where 


^ C3,fJ+-n  (3,f.)-n  (3,g>-n^2.82) 


as  before  one  defines  an  unperturbed  Green's  function  defined 
by 

F  fi,W  i)  = 

£'2 '3 '2'*'  -  ~  ~  ~  ~  ~  (2.83) 

<&(*»-?*;$(  3-3; 


*v 


Applying  the  same  analysis  as  described  in  section  IV,  one 
sees  that  G  9  can  be  approximated  as  t 


0fo  *  -  L  b  (P|-  ?l)  &  (fz~fl)  &  (2.84) 


-  32  - 


Therefore*  the  operator  Vc  for  this  case  can  be  approximated 
as 


Vc  •  <  V  Gv>  V  > 

b  $>(**-%)  Stf 3 -fj)  i  (■¥+-%) 

[n  ( j.fo + a  ( i  Jo -  a  ( ijj)-vGty <2-e5> 
df3df4cl| 


Upon  application  of  the  Markov  approximation*  one  obtains 


[2  A  (ho)  -A  (ht,-ts)-A(i,trb) 

-A  (  3.  —  A  (3,  fa'iO'M  ( (2-86) 

t  a  a,  <Ws)J 

By  (2. 28),’ one  then  obtains  the  equation  for  the  fourth  moment 
M4  «<  )  #<z,P,  )$>*(zff9)\j>%,U)  as. 

^  **  V  «V 


a  O 

(2.87) 
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vhere 

** 

=  2A(S.o)~  ACZJrf,)-  A(ZJrf*) 

-A  (l.  ffb)  -  A  (}■  fyf>)  t  A  <2*“’ 
+  A($,  U-h) 

V  ^ 

and 


A  (§.?*;=  2/c 


■J —CO 


(2.89) 
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Chapter  3  SOLUTIONS  FOR  THE  MUTUAL  COHERENCE 
FUNCTION 


I.  •♦Slow'*  and  "Fast"  Quasiparticle  Distribution  Functions 


Equation  (2.78)  is  a  single  mode  type  of  wave  problem 
dependent  on  the  statistics  of  the  random  medium.  In  order  to 
obtain  the  solutions  for  the  mutual  coherence  function*  we 
shall  define  the  following  Fourier  transform  with  respect  to 
the  "fast"  variable  £  *  xt  -  x%  ,  and  the  "slow”  variable 
x  »  xa  +  xj/2.  They  are  so  called  "quasiparticle  distribution 
function"  or  "Wigner  distribution  function"  described  independ¬ 
ently  by  Marcuvitz  (74)  and  Bremmer  (34).  One  introduces  the 
double  Fourier  representation. 


it  Mi -*»*»> 


ci  &t  <J 
(2K)* 


(3-1 ) 


A  A 


A  A 


where  ♦ia  $  (k^  ,  z)  and  V  t  (k*  ,  z)  are,  respective¬ 
ly,  the  Fourier  amplitudes  of  4^ 3  *P  (xi  »  z)  and  4»a  =*j*(xa  ,  z). 
Defining  the  "slow"  (  x  or  k  )  and  "fast"  (  fc  or  k  )  spatial 
and  wave  number  variables  via 


X  s 


*1  +  xa 

2 


t  9  *1  ~  X2 


K  a  “  ft* 


44  +  *t 
2 


(3.2) 
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whence 

~  $1*2  =  +  ** 


and 

=  jm* 


(3.3) 


One  rewrites  equation  (3.1 )  as 


2.JI  (3.4) 


where 


F  t  = 


/>  a  iKx  , 

<4,it>e  Ak 


2  71 


(3.5) 


is  called  the  "fast"  quasiparticle  distribution  function. 
Alternatively,  one  writes  equation  (3.1)  as 


<>*!,$.*>  =  j  S  <*'«.  5;  e  i  -&L 


(3.6) 


where 


<♦*?«*>  e*U-iA 


(3.7) 


and  termed  the  "slow"  quasiparticle  distribution  function 
also  called  the  ambiguity  function. 

The  reality  property  F(k,x,z)  *  F*(k,x,z)  of  the  fast 
quasiparticle  d.f.  follows  readily  from  eq. (3,5)$  however, 
positive  definiteness  is  not  in  general  assured.  S(K,f  ,}r) 
is  only  real  when  (  X )  is  an  even  function  which  can  be 
proved  using  eq.  (3.6).  S(K  , E  , £  ),  just  as  F(k,x,z)»  can  be 
negative  for  some  value  of  K  and  £  .  However,  it  should  be 

noticed  that  intergrals  of  f  F  ,  V  S  (*'$.$>  rrr 
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is  always  positive. 

It  is  interesting  to  observe  the  following  properties  of 
F(k,x,z)  and  S(  K  ,  E  ,  |  )« 

(a)  one  notes  that  from  eq . ( 3 . 4 )  and  (3.6) 

J  J  2ft  (3.8) 

which  suggests  that  S ( K  , E  »  J  )  can  be  treated  as  a  charac¬ 
teristic  function  of  F(k,x,z),  in  k,x  space,  and  vice 
versa.  The  use  of  S(  K  ,  E  ,  J  )  or  F(k,x,z)  mainly  dependent 
on  the  nature  of  the  problem  to  be  solved. 

(b)  for  E  *  0  and  K  «0  equations  (3.4)  and  (3.5)  yield 


JiPM  !•'  UMSJl.  -Hi'IWUlMl  !!i 
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=  [p (*-*' d 


(3.9b) 


whence  on  infers  that 


1  $  (*>3J|*>  J  *  =  >  2K 


(3.10) 


which  suggests  the  identification  of  F(k,x,z)  as  a  number 
density  of  "fast"  quasiparticles  in  k,x  phase  space. 

A 

(c)  Similary,  for  x  *  0  and  k  =  0  and  if>  (x)  is  even,  equations 
(3.6),  and  (3.7),  give 


(4,3)I2>=  J 


Sck.s,jj 


27t 


(3.11a) 


<l<£ 


S  IKS.})  <il 


(3.11b) 


Which  also  suggests  the  identification  of  S( «  , $  , ^  )  as  a 
number  density  of  "slow"  quasiparticle  in  «  ,  £  space, 
also, 


*r<%  »>J 
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=  jst«.t.S)4ia 


(d)  One  notes  that 


J 


=  S  (0,0-3)  s  N* (3) 


(3.13a) 


J  S  Ch.E.S 38  F  C°'0'  3)  5  (3.i3b) 


where  Nf(z)  and  Ns(z)  denote  the  total  number  of  "fast"  and 
"slow”  quasiparticles  respectively.  Both  Nf(z)  and  Ns(z) 
can  be  obtained  from  information  about  P  or  S. 

(e)  A  coarser  description  of  the  dynamics  of  a  quasiparticle 
system  views  the  overall  system  as  a  single  Macro-particle. 
It  is  of  interest  to  find  the  average  position*  momentum* 
size*  etc.,  of  the  overall  system.  The  average  coordinat¬ 
es  of  quasiparticles  are  defined*  in  terms  of  the  distribu¬ 
tion  functions  F  in  kfx  space*  by 


XCi) 


dfe_j£ 

an 


(3.14a) 


Similarly*  the  higher  moment  averages 


flK  vn  -  _  \ 

*  * 


*211  (3.16a) 


Kw£n-  = 


K^SCM,!)^ 

Ns(i)  '  '  2* 


(3.16b) 
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It  can  be  easily  proved  that 


m.  +  u 


3  (iE)n3(-*K)’ 


nS(K-J'5 


k=  l=o 


(3.17a) 


m+  ri 


K  2  3  C-ix)md  (ik)n  F 


(3.17b) 


I 


Therefore,  a  macro-description  of  the  overall  system  can 
be  found  if  either  F  or  S  is  known. 

II.  Exact  Solution  of  Mutual  Coherence  Function  for  the  Case 
H(  £  )  =  a  £*  (  cf.  Furutsu,  13) 

The  defining  equation  for  the  Mutual  coherence  function 
has  been  derived  in  section  (2.V).  Introducing  the  "fast"  and 
"slow"  coordinates 


f,  + 

X  =  — s _ x_ 

~  2 


l  =  f,  -  ft 


one  obtains 
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where  Ma(  x  ,  |  ,  }  )  *  ^  V  ( £i  >  3)  ^  (  j*  . 

We  now  consider  the  Fourier  transform  with  respect  tox.,  which 
is  exactly  the  "slow"  quasiparticle  d.f.  defined  in  section  I. 
The  resulting  differential  equation  for  S  (  K  »  6  »  5  )  becomes 

v  v 


+  H  U)  S  “  o 


(3.19) 


where 


Hce;  =  f  A  (°)  -  A  cij>J 


(3.20) 


Equation  (3.19)  has  the  form  of  a  collisional  "kinetic" 
equation  indicative  of  a  distribution  of  "slow  quasiparticle" 
in  K  ,  J  space  at  a  plane  z.  The  "Characteristic"  trajectory 
along  which  the  "particles"  move  is  defined  by 


=  O 


On  a  quasiparticle  trajectory,  one  observes  that 


44 =  - h  (B)  s  (K-t. 


(3.21) 


TON 
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or 

5  =  ScM^o^e  ° 

(3.23) 


This  implies  that  the  number  density  of  such  "slow  particles" 
per  unit  volume  of  the  phase  space  K  ,  £  decays  along  the 
trajectory  defined  by  equation  (3.23). 

For  waves  propagation  in  free  space,  the  defining  equa¬ 
tion  for  S  becomes 


S  (K,  to),  l)  =  S  (K (»JW,o) 

(3.24) 


which  provides  the  constancy  of  S  on  a  quasiparticle  trajectory. 

v 

To  illustrate  the  applicability  of  the  slow  q,p  view  we 
shall  consider  wave  propagation  through  a  random  medium  charac¬ 
terized  by  H( J )  =  a  5* •  As  an  initial  condition,  we  choose 
for  analytical  purposes  a  planar  beam  of  the  simple  form 


f(x,o)=  e 


(3.25) 


The  slow  quasiparticle  distribution  function  at  z  =  0  follows 
from  the  transform  relation  (3,6)  as 


43 


S  (*.£<  o) 


■1 


-  -^T  -aKX 

e  lk  e  dx 


d  X 


_.£i  .jsl^ 

b  e  4b‘  e  4 


(3.26) 


At  a  subsequent  distance  z  the  S  function  follows  from  (3.21) 
and  ( 3 . 23 )  as  * 

-<x)s  £*(3';<J 5' 

S  (.*,1.1)  =  S  Ck.  U«KO)e  ° 

-af  [S(°)  +  k}'J  J| 

=  s  (k.£c»),o^  e  J° 


-  5  (X.  8-k}'  °)  6  Jo 


-af  SC£-K(3-3';fd5' 

J  O 


2  I  2 


C£-Ksr  .  _o 

=  ji?be  4b*  e  4  * 

-$a[3£,*-3«K»x+  K'5^ 

.wie-"***'**'”  .. 
.  -  e*(-fe+a2)  *  ***£&♦  <ii) 


(3.27) 


For  the  beam  wave  case,  the  wave  at  the  aperture  z=0  can 
be  assumed  to  have  a  Gaussian  amplitude  distribution  with  the 
beam  size  W0  and  radius  of  curvature  R0  for  the  wave  front. 
Thus,  at  z=0  the  wave  is  given  by 


-4A.(«<r+  HOX* 

e  z 


(3.28) 


where 


(3.29) 


Introducing  the  normalized  coordinates 


one  rewrites  (3.28)  as 


-1 

f  (x.o)  =  e 


(3.30) 


Substituting  into  (3.6),  one  obtains  (for  simplicity,  we  will 
drop  the  bars) 
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sou.*;- Je  2 


=  Je 


-iK*  . 

e  d* 

-|t(r(;zV4‘J"'‘^x^  -*xx 

-  dx 


~±drS 


K ±1*1?  , _ 

^  J5 


(3.31) 


Using  equation  (3.23),  one  finds 


r=“  -  4o{r(£-K3J 

S(x.t.3)=  Pf-  e 


(3.32) 


note  that  (3.32)  can  be  reduced  to  (3.27)  by  setting  ofr= 
and  ot;,  =  0. 


(3.6) 


The  Mutual  coherence  function  at  z  can  be  calculated  via 

i,  <».!.»>•  j  5  («.«.»  e*K,-4r 


-Je 


-/lKl-BSN-C$K+iKX 


(3.33) 


v  v;,t. 

w'  S.  i .  'A'lA 
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where 


A  =  ( l  trfiY-aO;  +  3a^3 

4  <*r  ^ 


6  = 


_  i 


r  +  4-4^  +  °-l 


4  o< 


(3.34) 


C  »  4-r1  +  * 

2  2  «r 


whence  on  integration  of  (3.33),  one  obtains 


M.(*,£.  5)  * 


-  J*B  (ix+fc)1 


<2  n^A 


(3.35) 


where  A,B  and  C  are  defined  by  (3.34). 


Equation  (3.35)  is  the  exact  solution  for  the  Mutual 
coherence  function  at  a  distance  z.  It  is  of  interest  to  con¬ 
sider  as  a  function  of  Z  the  beam  size  ratio,  axial  intensity 
ratio,  and  the  normalized  MCF  at  x  =  0  for  the  following  cases t 
(case  a)  For  collimated  beams  (of**  0)  propagating  in  free 
space,  one  obtains  from  equation  (3.34) 

A  =  +  -4— 


4«<t* 


B  =  4-*r 


c  =  Y 


(3.36) 
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Then  the  beam  size  ratio  (E=0)  is  defined  by 


W.*(J>  _ 

4A 

Wc(o)  ~ 

=  i  +  3  V 


(3.37) 


increases  with  z. 

The  axial  intensity  ratio,  defined  as  the  intensity 
of  the  wave  at  x  =  £  =  0,  becomes 


I  (0.1)  _  .  1 

I  (0,0)  "  Ma  (.0,0,0)'  (1+aVJ*  „ 

( 3 • 38  ) 


and  decreases  with  z. 

The  normalized  Mutual  Coherence  function  (NMCF)  is  an 
important  function  for  the  measurement  of  coherence. 
The  integral  of  NMCF,  f  d  £  ,  gives  the 

value  of  correlation  scale.  At  x  *  0,  it  is  give  byi 


Mt(o.  t.j)  a 


°<rS*  /  '  1  \ 

4  'it  3Vr*  / 


(3.39) 


(case  b)  For  a  focused  beam  (*(^>0  )  and  a  divergent  beam 
Ua  <  0  )  propagating  in  free  space,  one  has 


mvf'  ■■r+t  s  WfPV* 
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6  *  4-^r  t 


(3.40) 


=  -*£■  -  4-  4*-  ( i  -  ^  3 ) 


2  ot 


The  beam  size  ratio  is  then 


W.  c  i) 

w.V; 


. 


,  U-^if 


*j.rix+ 


i 


(3.41) 


Whence  for  the  focused  beam,  a  minimum  beam  size 
occurs  at  a  distance  z  =  z '  such v that 


WoCy 

W }(*) 


=  0 


(3.42) 


Applying  equation  (3.41),  one  obtains 


}  * 


^  f 


(3.43) 


and  thus 


T’W 
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w.x<sj 

Wo  w 


_  (otWrf  (.  clS  \ 1 


(3.44) 


Using  equation  (3.33)  and  (3.38),  we  have  for  the 
axial  intensity  ratio 


IW 

1  (o) 


1 


ufr+u-'usf  (3*45) 


and  the  normalized  MCF  (x  *  0)  is  then 

C 

—  °(r  E*  f  j; _ "1 

*  Ur*3VCi-*i£j  (3*46) 

~  C 

We  note  that  the  correlation  width,  t  » 

is  minimum  when  z  ■  z*  ■  ■  . ^  , t  *  which  implies 

o*r  t 

that  the  correlation  length  of  a  focussed  beam  reduces 
to  a  minimum  at  the  focal  point  and  retains  its  co¬ 
herence  beyond  the  focal  distance. 

(case  c)  In  the  presence  of  turbulent  medium,  a>0,  the  beam 
size  ratio  becomes 


50  - 


W 

w«(») 


+  £l-o(r3  ^ 


(3.47) 


where  "a"  measures  the  strength  of  turbulence.  For 
the  collimated  beam,  dz  -  0,  the  beam  size  ratio 
reduces  to  ±  +  otraJz  +  ,  which  im¬ 

plies  a  broadening  of  beam  size  with  increasing 
distance.  If  «(* > 0 ,  indicative  of  focused  beams,  the 
minimum  values  of  (3.47)  will  occur  at 


3 


-  ~~  °^a,)  +  £  6>(r  T  df)  f  3  ^ 


4 


(3.48) 


If  rf;<0,  the  beam  focuses  at  negative  values  of  z, 
i.e.  it  diverges  with  increasing  z.  The  axial  inten¬ 
sity  ratio  is 


For  a  divergent  beam,  «CA’<0  ,  the  axial  intensity  ratio 
decays  with  increasing  turbulent  strength  or  increasing 
z.  For  a  focused  beam,  dz >  0  ,  I(z)/I(0)  will  reach  a 
minimum  value  at  z  *  z'  given  by  equation  (3.48). 


It.'J  'A'.' >; 


The  normalized  MCF  given  by  at  x  *  0  is  equal  to 


"  C 

=  e 


+  CLotri) 


Cl-UAf«r3J  <3-50) 


We  observe  that  for  collimated  beams  propagating  in  a  random 
medium,  wherein  3  0  and  a>0,  the  wave  will  loss  its  cohe¬ 
rence  (the  correlation  length  reduces)  with  increasing  z. 

For  *i>0,  the  focused  beam  will  reduce  its  correlation  length 
before  the  focal  point.  After  passing  the  focal  distance,  it 
broadens  until  turbulent  effects  dominate  at  which  point  the 
beam  again  loses  its  coherence.  For  the  divergent  beam  case, 
the  wave  will  increase  its  correlation  length  and  then  lose 
its  coherence  when  the  turbulent  effect  dominates. 

To  illustrate  the  evolution  of  Mutual  Coherence  Function 
plots  of  M(X  ,  £  ,  £  )  in  the  x  ,  £  space,  for  the  easels  1, 
a  ■  0  representation  of  a  focused  beam  propagating  in  free 
space,  are  shown  in  Fig.  3. la-3. lb  at  distance  z  3  0.5,  1,  1.5 
2.  For  a  initial  MCF  depicted  in  Fig.  1,  the  focusing  effect 
at  z  3  0.5  is  apparent. 

Figs.  2a-2b  display  the  Mutual  Coherence  Function  for  a 
focused  beam  in  a  random  medium  with  a  3  kg  3  1,  •***•(.* 3  1,  the 
beam  assumes  a  minimum  width  in  x  at  z*  according  to  equation 
(3.48)  and  then  spreads  in  the  x  direction.  Because  of  the 
turbulent  effect,  the  MCF  decorrelates  as  the  wave  is  propaga¬ 
ted  and  the  width  in  £  space  is  reduced.  The  amplitude  of 


Fig.  la. Amplitude  of  MCF  at  z  =  0.5;  a  focused 
beam  propagating  in  free  space;  ^  =  1. 
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Fig.  lc.  Amplitude  of  MCF  at  z  -  1.5;  a  focused 
beam  propagating  in  free  space;  «lT-  ct*  =  1. 
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Fig.  Id.  Amplitude  of  MCF  at  z  »  2;  a  focused 
beam  propagating  in  free  space;  el*.  !• 
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the  MCF  decays  with  increasing  z.  Normalized  MCFs  calculated 
by  use  of  equations  (3.38)  (3. 45) -(3. 50)  are  depicted  in  Figs. 
3-9.  For  the  case  of  a  collimated  beam  in  free  space,  the 
normalized  MCF  is  broadened  with  increasing  z  (Fig.3).  For 
the  focused  beam,  the  correlation  length  reduces  to  a  minimum 
value  at  z  =  0.5.  Beyond  the  focal  point,  the  be  aim  starts  to 
spread  and  increases  its  correlation  length  (Fig. 4).  Fig.  5 
displays  the  normalized  MCF  for  a  collimated  beam  in  a  medium 
with  the  turbulent  strength  a  =  1}  the  correlation  length 
decreases  as  z  increases.  The  same  case  for  a  focused  beam 
is  depicted  in  Fig.  6.  Due  to  focusing,  the  width  of  norma¬ 
lized  MCF  becomes  narrower  for  small  z,  it  then  broadens  for 
some  distance  because  of  the  diffraction  effect,  finally, 
turbulent  effects  reduces  its  coherence  length.  Fig.  7  and  8 
shows  the  MCF  different  values  of  "a”  for  a  collimated  beam  and 
a  focused  beam  respectively.  In  both  cases,  the  MCF  become 
more  decorrelated  as  the  turbulence  becomes  stronger.  The 
effects  of  turbulence  for  different  initial  beam  sizes  are 
shown  in  Fig.  9a-9b  for  collimated  beam  and  focused  beam.  The 
correlation  length  is  larger  for  a  wider  beam  size  ,  as  one 
expects .  i 

III.  Numerical  Evaluation  of  MCF  for  a  Complicated  Power 
Spectrum 

To  obtain  MCF  for  an  arbitrary  power  spectrum,  one  can 
apply  the  same  procedure  as  described  in  section  II.  However, 
because  of  the  complicate  integration  involved  in  evaluating 
the  function  H(g),  it  is  not  realistic  to  solve  the  problem 
analytically.  A  numerical  procedure  via  the  quasiparticle 
distribution  function  approach  will  be  used  to  find  the  d is ire 
information. 

It  is  necessary  to  develope  a  numerical  scheme  which  cal¬ 
culates  the  slow  quasiparticle  distribution  function  S(K  ,£  ,5) 
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Fig.  2a.  Amplitude  of  MCF  at  z  =  0.5;  a  focused 
beam  propagating  in  a  random  medium  with  H(£)= 
a  £ 1  ;  a  =  1,  o<2  =  oC  =*  1. 
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Fig.  2b.  Amplitude  of  MCF  at  z  =  lj  a  focused 
beam  propagating  in  a  random  medium  with  H(£) 
a  £*>  a  *  1,  o ir  -  =  1. 
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Fig.  5.  Normalized  MCF  vs.  fc  for  different  values 
of  zf  a  collimated  beam  propagating  in  a  random 
medium  with  H(fc)  =  a  a  =  1. 
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Fig.  6.  Normalized  MCF  vs.  £  for  different  values 
of  zj  a  focused  beam  propagating  in  a  random  medium 
with  H(|)  »  at1,  a  •  1. 
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Fig.  7.  Normalized  MCF  vs.fc  for  different  values 
of  a)  a  collimated  beam  propagating  in  a  random  medium 
with  H ( £  )  =  a  £ 1 . 
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of  a |  a  focused  beam  propagating  in  a  random  medium 
with  H(  &  )■  a  . 


Fig.  9a.  Normalized  MCF  vs.  £  for  different  values 
of  beam  width  b;  a  collimated  beam  propagating  in 
a  random  medium  with  H(  £  )  =  a  £*  .  a  a  1. 
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Fig.  9b.  Normalized  MCF  vs.  E  for  different  values 
of  beam  width  b|  a  focused  beam  propagating  in  a 

random  medium  with  H(  t  )  *  at1',  a  «  1. 
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for  a  given  S  (  K  »  &  ,  6  ) .  The  corresponding  results  for  inten¬ 
sity,  beam  size  ratio  and  normalized  MCF  can  then  be  obtained 
via  equations  (3.6)  and  (3.11).  We  shall  first  examine  the 
case  of  a  collimated  beam  propagating  in  free  space  to  check 
the  accuracy  of  the  algorithm  being  used.  The  z  =  0  quasipar¬ 
ticle  distribution  function  which  follows  from  equation  (3.31), 

1  e  ©  e 


1  ,  _  (K +«**£)* 

s  e'T  rE  e"  'Mr 

(3.51) 


is  depicted  in  Fig.  10  for  the  case  *r=  1,  =0.  The  quasi¬ 

particles  are  evenly  distributed  in  the  intervals  -7  <  £  <  +7, 
-7<K<  +7.  At  a  subsequent  distance  z  the  distribution  func¬ 
tion  follows  from  (3.24)  and  (3.21)  as 


So*,*,?*;  =  S  <*/ ° ) 

(3.52) 


i.e.  the  quasiparticle  density  is  redistributed  when  evolutes 
in  ,  as  the  quasiparticle  spread  in  £  the  probability  dis¬ 
tribution  function  (d.f.)  with  coordinates  S(  K,  0),  at 

z  *  0  "moves”  to  the  new  location  S(K  ,  £  ,  J'1)  at  z  =  z*.  The 
redistribution  phenomena  enable  us  to  interpolate  every  point 
in  K  ,  £  space  from  the  initial  d.f.  to  the  final  d.f.  at  z*»z* 
numerically.  Figs. 11-12  display  the  S  function  at  z  ■  0.75, 

1.5  respectively.  The  quasiparticles  with  higher  K  "move" 
faster  than  those  with  smaller  K  .  The  resulting  d.f.  thus 
spreads  along  £  directions.  The  MCF  at  x  *  0  as  a  function  of 
£  for  distances  z  corresponding  to  the  associated  S (*,£,£  ) 
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Fig.  10.  Initial  slow  quasiparticle  distribation 
function  for  collimated  beam  with  e<r  8  li  z  »  0. 
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Fig.  11.  S  function  in  K  ,  £  space t  z  ■  0.75 
a  collimated  beam  propagating  in  free  space. 
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Fig.  12.  S  function  in  K  ,  £  space;  z  *  1.5, 
a  collimated  beam  propagating  in  free  space. 
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Fig.  13.  MCF  vs.  2  at  x  ■  0;  a  collimated  beam 
propagating  in  free  space. 
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Fig.  19.  Intensity  function  I(x,z)  vs.  x 
a  focused  beam  propagating  in  free  space 
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are  displayed  in  Fig.  13.  The  spreading  effect  increases  the 
coherence  length  as  z  increases.  Fig.  14  depicts  the  initial 
condition  described  by  equation  (3.51)  with  elr  =  =  1,  a  =  0 

i.e.  focused  beam  in  free  space  at  z<0.5,  the  field  reduce 
its  coherence  length  because  of  the  focusing  effect,  it  then 
spreads  as  z  =  0.5,  as  described  in  section  II.  The  S  func¬ 
tion  at  z  *  0.75,  1.5  are  shown  in  Fig.  15,  16.  Fig.  17  dis¬ 
plays  the  MCF  as  a  function  of  £  and  z.  The  above  exact  re¬ 
sults  have  been  used  to  check  the  accuracy  of  the  algorithm 
and  found  to  be  very  good. 

The  Fourier  transform  of  S( H  , £  ,  ^ )  gives  the  MCF 
M(x  ,  £  ,  ).  Setting  £  =  0,  one  obtains  the  intensity, 

I(X  ,  *  )  =  M(X,  0,  5-  )  as  a  function  of  x  and  z.  Figs.  18 
and  19  depict  the  intensity  function  for  collimated  and  focused 
beams  propagating  in  free  space.  In  the  first  case  the  inten¬ 
sities  spreads  in  the  X  direction  as  z  evolves.  The  focusing 
effect  is  observed  in  the  second  case  and  is  shown  in  Fig. 19. 

For  a  wave  propagating  thru  a  random  medium,  characteriz¬ 
ed  by  a  power  spectrum  £  (A  )$  the  evolution  of  the  S  function 
can  be  obtained  from  the  relations 

5  (k.e.I;  =  S  (K .£<•>. o;e  0 

(3.53) 


where,  for  the  one  dimensional  case,  H(£  )  can  be  represented 
as 


n  co 

HU)*  \m  Cl  -  COS  ill  )  $  (At)  dfc 


(3.54) 
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Along  the  quasiparticle  trajectory,  one  can  rewrite  (3,53)  as 


_  j5H 


5  ck,£,5  j)  =  S  ck,  s-Kj,o)e 


(3.55) 


From  knowledge  of  the  power  spectrum  and  the  initial 
distribution  function,  one  can  interpolate  the  S  function 
along  the  trajectory  in  accordance  with  (3.55).  That  is,  the 
S  function  spreads  in  £  space  and  decays  along  the  trajec¬ 
tory  with  a  structure  function  H( £  )  dependence. 

One  finds  equation  (3.55)  is  different  from  eq.  (3.52) 
by  a  collisional  (decaying)  term  exp 

It  is  thus  important  to  observe  the  behavior  of  the  H  func¬ 
tion  for  an  arbitrary  power  spectrum.  We  shall  choose  three 
different  kinds  of  power  spectrum  and  study  the  resulting 
behavior  of  the  H  function! 

(a)  Gaussian  spectrum 

_  (»-*’> 


where  is  the  center  of  the  Gaussian  spectrum  and  b  is 
the  width 

(b)  Power  law  spectrum  with  spectral  index  P 


^  (6)  =  Ap  ( 1  +  ^  z  ) 


(3.57) 
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where  1/l©  is  the  wavenumber  associate  with  the 

outer  scale  L0  and  P  is  the  spectral  index. 

(c)  Compound  spectrum 

/  o  D>\* 

$  (V  =  ^  [(1  +  ^  + 

is  the  superposition  of  spectrum  (a)  and  (b). 

The  power  spectrum  and  associated  H  function  are  shown 
in  Figs.  20a  and  20b.  Curve  1  in  Figs.  20a,  b  represents  $  (*) 
and  H(£  )  for  spectrum  type  (b).  The  Gaussian  bump  power  spec¬ 
trum  and  corresponding  H  function  are  depicted  in  curves  2,3,4 
with  ik'  equal  to  1,2,3  respectively,  and  with  a  =  0.3,  b  *  0.2. 
One  observes  that  the  additive  Gaussian  bump  contributes  an 
oscilltory  structure  to  the  power  spectrum.  For  instance,  •&’  =3 
imposes  an  oscilltory  structure  with  three  peaks  located  at 
£  =  Jt/3,  7L  ,  5  It/3. 

The  oscillatory  behavior  implies  a  higher  randomness  in 
the  medium.  As  a  result,  the  wave  is  more  decorrelated  by  the 
random  field  and  a  stronger  scattering  effect  is  encountered. 
Fig.  21a  displays  the  Amplitude  dependence  of  the  Gaussian 
bump  for  a  =  0.3,  0.5  and  0 , 7  from  lower  curve  to  higher  curve . 
The  superposed  behavior  of  spectrum  (a)  and  (b)  is  shown  in  Fig. 
21b  for  the  case^'-l,  a  =  0.3,  b  -  0.2. 

It  is  of  interest  to  evaluate  numerically  the  S  function 
at  a  distance  J  using  equation  (3.55).  Once  the  S  function 
is  obtained,  one  can  easily  find  the  intensity  function  I(x,z) 
and  correlation  function  M(0, £  ,  J  )  at  x  =  0  from 


a  e 


2  b* 


(3.58) 


OIST*  2.13281  WIDTH*  1.08008  Al«  1.80000  GAMMA*  1 . 00000 
GCTR*  3.00000  GUID*  2.00000E-01  GAMP*  3.00000E-01  INDEX*  2.00000 
POWER  SPECTRUM  WITH  A  GAUSSIAN  BUMP 


DIST*  2.13281  WIDTH*  1.00000 

GCTR*  3.00000  GW ID*  2.00000E-01 

POWER  SPECTRUM  WITH  A  GAUSS  I AH  BUMP 


AI*  1.00000  GAMMA*  1.00000 
GAMP*  3.00800E-01  INOEX*  2.08000 


DIST*  2.13281  WIDTH*  1.88098  AI*  1.88808  GANNA*  1.89888 
GCTR*  2.88988  GH10*  2.08088E-81  GANP*  7.08880E-81  INDEX*  2.88889 


8.8  1.25  2.38  3.73  3.88 


a 

Fig.  21a.  H(  £  )  vs.  J  ;  spectrum  (c)  with  a  =  0.3 
(curve  1),  a  *  0.5  (curve  2  ),  a  =  0.7  (curve  3). 
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Fig.  21b.  H(  t  )  vs.  fc  j  curve  1  for  spectrum  (a), 
curve  2  for  spectrum  (b),  curve  3  for  spectrum  (c) 
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The  correlation  scale  &  of  the  field  is  defined  by 
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and  hence  the  knowledge  of  M  (0,  £  ,  z)  permits  one  to  eva¬ 
luate  the  correlation  scale  1.  We  shall  first  present  some 
numerical  evaluations  of  the  S  function  by  means  of  eq.  (3.55). 
The  significance  of  intensity  and  correlation  length  for  both 
focused  and  collimated  beams  will  be  discussed  later. 

In  the , following  numerical  evaluations,  we  have  used  the 
normalized  quantities i  z  -♦  zA«L,'  ,  x -*  x/L»  ,  V  =%kf  L*A(0). 
where  ka  is  the  wavenumber  of  the  incident  field  and  L«  is 
the  outer  scale  of  the  random  field.  Figs.  22a  and  22b  depict 
the  S  function  at  the  boundary  z  *  0  for  the  collimated  beam 
and  focused  beam,  respectively.  Figs.  23a-d  depict  the  S  func- 
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Fig.  22a.  S  function  Bor  a  collimated  beam  with 
of»-  *  1}  z  =  0. 
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Fig.  23a.  S  function  for  a  focused  beam  propa¬ 
gates  thru  a  random  medium  with  a  power  law  spectrum 
p  =  4,  =  10,  otr  =  <**'  =  1  >  z  =  0.4375. 
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tion  at  z  =  0.4375,  0.875,  1.3125,  1.75  for  focused  beam  pro¬ 
pagating  thru  a  random  medium  with  power  law  spectrum  defined 
by  the  spectral  index  p  =  4,  f  =  10,  *  1.  One  observes  that 

the  quasiparticle  d.f.  spreads  in  the  E  direction  and  decays 
according  to  the  H  function  as  it  redistributes.  If  a  Gaussian 
bump  is  imposed  with  its  center  located  at  k'  -  3(GCTR),  width 
b  *  0.2(GWID),  amplitude  a  =  0.3  (GAMP),  the  S  function  is  dis¬ 
torted  and  sharpened  compared  to  the  case  of  the  simple  power 
law  spectrum.  These  effects  are  displayed  in  Figs.  24a-24d. 

As  discussed  earlier,  the  location  of  the  Gaussian  bump  deter¬ 
mines  the  number  of  peaks  in  the  H  function.  In  order  to  illu¬ 
strate  the  effect  of  these  peaks  on  the  S  function,  we  shall 
examine  the  case  of  a  collimated  beam  propagating  thru  a  random 
medium  with  the  composite  spectrum  (c).  The  center  of  the 
bump  is  chosen  to  be  located  at  k*  *  1,  3,4,  respectively.  For 
a  10,  GWID  =  0.2,  GAMP  *  0.3,  p  =  4,  one  displays  the  nume¬ 
rical  solutions  in  Figs.  25a-25d,  figs.  26a~26d  and  figs.  27a- 
27d.  At  z  =  0.4375,  the  S  function  is  less  distorted  for  k'=l 
than  it  is  for  k*’  =  3  and  4.  At  z  =  1.75,  the  S  function  for 
k'  =  3  and  4  clearly  shows  the  effects  of  the  H  function  peaks. 
One  observes  that  the  presence  of  these  peaks  in  the  H  function 
greatly  reduces  the  correlation  length  in  the  E  direction. 

It  is  easier  and  clearer  to  discuss  the  independent 
effects  of  different  spectrum  structure  by  means  of  the  inten¬ 
sity  function  and  correlation  length  associated  with  the  propa¬ 
gating  waves.  Substituting  the  numerically  evaluated  S  func¬ 
tion  into  equations  (3. 59)- (3.60) ,  one  obtains  the  desired 
information  for  I(x,z)  and  £  .  Fig.  28a  depicts  collimated 
beams  propagating  thru  a  random  medium  with  H(  E  )■  VS*  . 

Curve  1  represents  the  axial  intensity  and  correlation  length 
for  the  case  of  wave  propagation  in  free  space.  Turbulent 
effects  for  T *  0.5  and  1  are  depicted  in  curves  2  and  3.  (In 
Figs.  28-31  we  use  the  same  curve  index  for  different  values  of 
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Fig.  24a.  S  function  for  a  focused  beam  propagates 
thru  a  random  medium  characterized  by  spectrum  (c)t 
k'=  3,  b  a  0.2,  a  =  0 . 3 ,  ok  »  =  1»  z  =  0.4375. 
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Fig.  24c.  2  =  1.3125. 
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Fig.  25a. S  function  for  a  collimated  beam  propagates  thru 
a  random  medium  characterized  by  spectrum  (c)jk’  »  1, 
a  *  0.3,  b  »  0.2,  p  »  4,  °f  -  10}  z  -  0.4375. 
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Fig. 26a. S  function  for  a  collimated  beam  propagates  thru 
a  random  medium  characterized  by  spectrum  (c);  k'  *  3, 
a  *  0.3,  b  a  0.2,  p  =  4,  -  10,  z  =  0.4375. 
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Fig.  27a.  S  function  for  a  collimated  beam  propagates 
thru  a  random  medium  characterized  bt  spectrum  (c); 

*  4,  a  *  0.3,  b  =  0.2,  p  *  4,  10,  z  *  0.4375. 
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y  ).  One  observes  that  the  axial  intensity  1(0, z)  decays  as 
the  wave  propagates  while  the  correlation  length  decreases  due 
to  turbulent  effects  until  about  z  =  0.3;  diffraction  effects 
then  become  important  and  coherence  is  partially  restored. 

After  further  propagating  turbulent  effects  again  dominate  and 
the  wave  is  decorrelated.  For  a  focused  beam  propagating  in 
free  space  wither*  1,  &L-1,  the  intensity  at  z  =  0.5  will 
reach  a  maximum  and  its  correlation  length  a  minimum  as  shown 
in  the  curve  1  of  Fig.  28b,  If  a  turbulent  medium  is  present, 
the  focal  distance  is  reduced  due  to  the  scattering  effect  im¬ 
posed  by  the  turbulent  medium.  Curves  2,3  of  Fig.  28b  display 
the  decorrelation  process  contributed  by  both  focusing  and 
turbulent  phenomena  from  z  =  0  to  z  s’  0.4.  At  z  =  1.2,  the 
turbulent  effect  dominates  and  reduces  the  correlation  length 
as  z  increases.  Figs.  29a, b  display  the  same  effect  for  the 
Gaussian  spectrum.  Comparing  these  curves  with  the  square  law 
spectrum,  one  observes  that  the  Gaussian  spectrum  has  a  small¬ 
er  effect  on  the  wave  statistics  than  the  square  law  spectrum. 

It  is  of  interest  to  consider  a  power  law  spectrum  for  the 
ionospheric  random  medium.  We  shall  choose  the  spectrum  type 
(b)  with  p  =  4  and  make  use  of  numerical  calculations  of  the 
S  function  to  evaluate  the  axial  intensity  function  I(o,z)  and 
correlation  length  1.  The  results  are  depicted  in  Figs.  30a, b. 
Fig.  30a  displays  the  case  of  collimated  beam  for  -  0.5  and  1. 
The  turbulent  effect  slightly  reduces  the  intensities  (Note  that 
this  is  a  weak  turbulence  case),  and  the  correlation  length  de¬ 
creases  compared  to  the  case  of  free  space  (curve  1).  For  the 
case  of  focused  beams,  depicted  in  Fig.  30b,  one  finds  that 
scattering  effects  decrease  the  focal  length;  the  wave  is  also 
diffracted  and  decorrelated  as  in  the  previous  case.  The  spec¬ 
trum  of  a  real  random  medium  is  not  a  smooth  function  in  general 
Experimental  results  show  that  the  power  spectrum  could  posses 
a  number  of  "bumps"  in  various  frequency  ranges.  We  shall  illus 
trate  the  effect  of  such  bumps  by  simply  adding  a  typical 


0r3T*  2  13251  WIDTH*  1.00000  01*  0.0  GAMMA*  00 

GCTR*  0  0  GW 10*  1.00000  GAMP*  1.00000  INDEX*  2.00000 


Fig.  28a.  -[Hi  1  vs.  z  dotted  line;  1(0, z)  vs.  z 
solid  line;  curve  li  free  space;  curve  2j*T’  =  0.5; 
curve  3i  Y  *  1;  collimated  beam  with  «*>- =  1. 

H(  l  )  =  r  E  . 

OIST*  2  13281  WIDTH*  1.00000  AI«  1  30000  GAMMA*  3  00000E-01 

GCTR*  0.0  GUID*  1.00000  GAMP*  l.:00000  INDEX-  2.00000 


■> 


Fig.  '28b.  focused  "beam  case;  o(w  =  o(:  *  1, 


D!ST>  2.13281 
GCTR*  1.1 
GAUSSIAN  SPECTRUM 


UIOTH-  1.08808 
CNIO-  1.80088 


AI-  M 
GAMP-  1.80808 


GANNA-  0.8 
INDEX*  2.08008 


. . 


5.00E-01 


2 

. 


0.0  3.33E-01  1.0?  1.S0  2.13  ‘ 

2 

Fig.  29a.JziL  1  vs.  z  dotted  line;  1(0, z)  vs.  z 
solid  line;  (1)  free  space;  (2)  r»  0.5j  (3)«y 
Gaussian  spectrum,  k’  =0,A=l,b=l. 


OIST*  2. 13281 
GCTR-  0.0 
GAUSSIAN  SPECTRUM 


WIDTH-  1.00000 
GU10-  1.00000 


AI-  1.00000 
CAMP-  1:00000 


GAMMA-  1.00000 
INDEX-  2.00000 


i, is l 


5.80E-01 


3.33E-81  1.17 


Fig.  29b.  focused  beam  case;  oLr  *  -  l. 
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Gaussian  bump  on  the  power  law  spectrum,  i.e.  a  compound  spec¬ 
trum  as  described  in  eq.  (3.58).  Fig.  31a,  31b  depict  the 
effect  of  a  Gaussian  bump  on  the  axial  intensity  function  and 
the  correlation  length  for  collimated  and  focused  beams  res¬ 
pectively.  The  parameters  are  chosen  as  t  T  -  0.5  and  1. 

GCTR  ■  1,  GWID  =  0.2,  GAMP  =  0.3  and  p  =  4.  The  results  show 
that  the  addition  of  this  Gaussian  bump  causes  strong  scatter¬ 
ing  phenomena.  We  have  presented,  for  the  weak  fluctuation 
cases  V  =  0.5,  1,  the  effects  of  turbulence  in  the  medium  . 

As  a  comparison,  we  shall  now  display  in  Figs.  32-33  some  re¬ 
sults  for  the  strong  turbulence  case,  ~  5  and  10.  A  colli¬ 
mated  beam  propagates  thru  a  random  medium  with  a  composite 
spectrum  formed  by  a.  Gaussian  bump  located  at  GCTR  =  1,  GWID  = 
0.2  GAMP  =  0.3,  for  the  cases  "T  -  5  and  10.  Curves  2,3  in  Fig, 
32a  represent  the  intensity  and  correlation  length  for  -  5 
and  10;  they  shows  that  for  strong  turbulence  case  the  corre¬ 
lation  length  decreases  rapidly  and  a  long  tail  appears  at 
large  Z.  The  axial  intensity  drop  appreciably  and  approachs 
a  small  constant  value.  For  the  focused  beam  case,  the  focal 
length  is  greatly  reduced  for  large  values  of  *T  .  We  shall 
expect,  as  the  turbulent  strength  is  very  strong,  that  the 
focusing  effect  will  not  be  present  and  the  wave  will  decay 
as  z  increases.  These  results  are  depicted  in  Fig.  32b.  If 
a  Gaussian  bump  is  located  at  k’  *  3,  the  scattering  effect  is 
stronger  than  for  the  previous  case.  Figs.  33a  and  33b  dis¬ 
play  data  for  a  collimated  beam  and  focused  beam  propagating 
thru  a  random  medium  described  by  a  composite  spectrum  with 
k'  »  3.  One  observes  that  decay  is  much  greater  for  the  in¬ 
tensity  and  that  decorrelation  phenomena  are  more  evident 
than  for  the  case  k*=l  . 
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b  =  0.2,  p  =  4.  collimated  beam  case. 
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Fig.  32a.  Hit  1  vs.  z  dotted  line;  l(0,z)  vs.  z 
solid  line;  (1)  free  space;  (2)  Ts  5;  (3)  T  =  1 
Power  spectrum  with  a  Gaussian  peak,  k'  =  1,  a  =  0 
b  *  0.2,  p  =  4,  collimated  beam  case. 
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IV.  Conclusion  and  Discussion 


A  wavepacket  or  beam  can  be  viewed  as  a  finitely  extend¬ 
ed  moving  and  deforming  system  of  "Quasiparticles".  The  Four¬ 
ier  transform  of  the  Mutual  Coherence  Function  w.r.t.  either 
the  "fast"  or  "slow"  coordinate  is  defined  respectively  as  the 
"fast"  or  "slow"  quasiparticle  distribution  function  F  or  S. 

For  wave  propagation  in  a  random  medium,  it  is  usually  more 
convenient  to  utilize  the  "slow  quasiparticle"  distribution 
function  S.  The  knowledge  of  the  S  function  provides  the  de¬ 
sired  information  for  two  point  statistics  of  a  stochastic 
wave.  For  example,  the  correlation  length  lat  any  value  z, 
defined  in  equation  (3.61),  can  be  interpreted  as  the  total 
•  numer  of  q.p.  Ns(z)  divided  by  the  q.p.  density  at  E  =  0. 

Also,  the  axial  intensity  l(0,z)  may  be  represented  by  the 
slow  q.p.  density  at  £  =  0.  To  obtain  the  2  point  statisti¬ 
cs  at  any  distance  z,  one  has  only  to  observe  the  q.p.  flow 
in  K  ,  E  space  as  z  increases.  For  the  case  of  a  collimated 
beam  propagating  in  free  space,  one  considers,  at  z  =  0,  a 
system  of  quasiparticles  evenly  distributed  between  -7  <  £  , 

K  <  +7.  Qualitatively,  we  shall  use  the  "area"  of  a  circular 
"point"  to  represent  the  number  density  of  slow  q.p.;  thus  at 
z  -  0  the  q.p.  are  distributed  in  K  »  £  space  as  shown  in 
Fig.  34.  As  z  increases,  those  q.p.  with  higher  K  "move" 
faster  than  those  with  smaller  K  .  As  a  result,  the  q.p. 
density  at  £  =  0  decreases  as  z  increases  and  hence  the  axial 
intensity  1(0, z)  decreases.  Since  the  total  number  of  q.p  is 
conserved  in  free  space  propagation  ,  i.e.  Ns(z)  is  constant, 
the  correlation  length  increases  as  q.p.  move  along  the  trajec¬ 
tory. 


For  a  focused  beam  propagating  in  free  space,  the  q.p. 
at  z  =  0  are  distributed  in  K  »  £  space  with  a  maximum  along 
an  axis  deviating  from  the  £  =  0  axis  as  shown  in  Fig.  35. 
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One  observes  that  the  q.p,  density  reaches  a  maximum  at  a 
plane  z  -  Zf,  when  those  q.p.  with  larger  density  reach  the 
center  line  6=0.  i.e.  the  axial  intensity  is  maximum  at 
the  focal  distance  Zf.  The  conservation  of  the  total  num¬ 
ber  of  q.p.  in  free  space  then  implies  a  reduction  of  the 
correlation  length  from  z=0toz=Zf.  As  the  q.p.  flow 
moves  beyond  the  focal  point,  those  q.p.  with  lower  density 
move  to  the  center  and  the  axial  intensity  decreases  as  z 
increases;  Correspondingly,  this  implies  an  increase  in  the 
correlation  length. 

For  a  wave  propagating  in  a  random  medium  characteri¬ 
zed  by  the  structure  function  H(  6  ),  one  just  simply  leads 
to  a  decay  of  the  wave  field  in  the  K  ,  6.  space  as  depicted 
in  Fig.  36.  The  structure  function  acts  as  an  "absorber" 
which  absorbs  part  of  the  q.p.  according  to  its  strength. 

For  z>0.  Fig.  37  depicts  for  a  collimated  beam  at  z  =  0, 
how  the  q.p.  move  in  a  random  medium  characterized  by  H( 6  ). 
As  z  increases  those  q.p.  with  smaller  density  move  toward 
the  center  while  some  of  them  are  absorbed  by  the  structure 
function.  The  resluting  q.p.  density  at  6 =  0  is  decreased 
and  is  smaller  than  in  the  case  of  free  space.  One  notes 
that  the  total  number  of  q.p.  is  not  conserved  in  this  case. 
The  ratio  of  Ns(z)  and  l(0,z),  which  determines  the  correla¬ 
tion  length  1,  will  decrease  as  z  increases  for  the  strong 
turbulent  case;  however,  in  the  weak  turbulence  case,  it  will 
first  decrease  and  then  increase  as  z  increases  at  least  for 
a  short  distance  (Note  that  the  q.p.  at  K  =  0  will  never 
move).  For  the  case  of  a  focused  beam  propagating  in  a  ran¬ 
dom  medium,  one  depicts  the  q.p.  flow  thru  a  random  medium 
characterized  by  H(  £  )  as  in  Fig.  38. 

Because  of  turbulence  in  the  medium,  the  q.p.  vith 
larger  density  are  absorbed  before  reaching  the  center  £  =  0. 
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As  a  result,  the  q.p.  density  at  £  *  0  creates  a  maximum  as  a 
result  of  q.p.  which  initially  have  less  density.  Thus,  the 
focal  length  (zf)  is  smaller  compared  to  the  free  space  case. 
The  correlation  length  in  this  case  is  reduced  from  z  =  0  to 
z  =  Zf.  Beyond  the  focal  point,  the  q.p.  density,  at  £  =  0, 
is  reduced  and  the  ratio  of  Ns(z)  and  l(0,z),  which  yields  1, 
will  again  decrease  for  the  strong  turbulent  case  or  decrease 
and  then  increase,  as  z  increases  at  least  for  a  short  dis¬ 
tance,  for  the  weak  turbulent  case. 

The  qualitative  discussion  presented  above  is  suitable 
for  an  arbitrary  random  field.  The  form  of  the  structure 
function  H(  £ )  will  determine  the  q.p.  density  profile.  To 
illustrate  this  effect,  we  will  present  some  numerical  re¬ 
sults  for  a  wave  propagating  in  a  random  medium  characterized 
by  the  spectrum  (C).  Fig.  39a  depicts  the  q.p.  density  func¬ 
tion  at  x  =  0  in  £  -space,  i.e.  M(0,  £  ,  z),  for  the  case  of 
collimated  beams  propagating  in  a  turbulent  medium  with  spec¬ 
trum  (C).  We  choose  the  parameters  GCTR  =3,  GWID  =  0.2, 

GAMP  =  0.3  and  p  =  4.  One  notes  that  the  axial  intensity 
decreases  as  z  increases  and  the  deformation  of  M(0,  £  ,  z) 
reflects  the  behavior  of  the  decay  field.  The  normalized 
correlation  function,  | j  is  shown  as  a  function  of  z 

in  Fig.  39b.  The  area  under  each  curve  represents  the  corres¬ 
ponding  correlation  length  1.  For  T  -  10,  which  represents 
a  strong  turbulent  case,  one  sees  that  the  correlation  width 
decreases  as  z  increases. 

The  corresponding  phenomena  for  the  focused  beam  case 
is  displayed  in  Fig.  40a  and  40b.  One  notes  the  focusing  and 
defocus ing  effects  in  Fig.  40a  and  the  associated  curves  for 
the  normalized  MCF  in  Fig.  40b. 

The  choice  of  an  symmetric  function  for  the  input  field 
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Fig.  37.  S  function  of  a 
collimated  beam  moves  in  a 
random  medium. 


Fig.  38.  S  function  of  a 
focused  beam  moves  in  a  random 
medium. 


OIST*  1.75000 
GCTR*  3.00000 


MOTH*  1.00000 


AI*  0.0 


GAMMA*  10  0000 


GMD*  2.00000E-01  GAMP*  3.03000E-01  INDEX*  2.00000 


POWER  SPECTRUM  UITH  A  GAUSSIAN  8UMP 
■X 


Fig.  39a.  q.p.  density  function,  at  x  =  0  ,  in  £ 
space;  plotted  as  a  function  of  £  and  z;  collimated 
beam  propagating  thru  a-  random  medium  with  spectrum  (c) 
k'  ■  3,  a  =  0.3,  b  =  0.2,  p  =  4,  -r  =  10. 
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Fig.  39b.  Normalized  q.p.  density  function,  M 
vs.  £  and  z;  the  area  under  each  curve  represents 
the  corresponding  correlation  length  1. 
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Fig.  40a.  q.p.  density  function,  at  x  =  0,  in  t  space; 
focused  beam  propagating  thru  a  random  medium  with 
spectrum  (c),  k*  =  3,  a  =  0.3,  b  =  0.2,  p  =  4,  *y*  =  10. 
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leads  to  a  real  S  function.  It  provides  a  simpler  approach  to  l 

obtaining  2  point  statistics  especially  for  numerical  evalua-  [ 

tions.  The  exact  solution  of  the  Mutual  Coherence  Function  I 

for  the  case  H(  £  )  =  a  is  solved  and  presented  in  section 

II.  For  an  arbitrary  spectrum,  a  numerical  scheme  has  been 
used  to  evaluate  the  S  function  as  shown  previously  in  section 

III. 


For  simplicity  the  above  analysis  has  been  limited  to 
the  two  (z,x)  dimensional  case.  The  extension  to  three  dimen¬ 
sions  is  straightforward  but  requires  more  computation.  For 
a  random  medium  which  is  inhomogeneous  in  z,  one  simply  adds 
a  z  dependence  to  the  H  function  in  equation  (3.55)  and  all 
of  the  calculational  procedures  will  then  be  ^unchanged . 
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Chapter  4  THE  FOURTH  ORDER  MOMENT  AND  INTENSITY 
FUNCTION 


I.  Formulation  of  the  Fourth  Moment  Equationt 


The  parabolic  equation  for  the  4th  moment 

m,  ( r. ,  r. .  ?•')  =  <■?  w  v  <m.)  y'yt.') 

^  V  v  / 


(4.1) 


may  be  derived  for  the  case  of  statistically  homogeneous  fluc¬ 
tuation  as  follows s  (cf.  Tatarskii,  1;  Ishimaru,  25) 
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where 
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and 

H  i3,f.)  =  A  (  3-0  -  A  (i-f) 


(4.4) 


Introducing  the  new  variables! 
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n  = 


then  Vp*  +  Vf4  -  V^<  -  =  2(  Vp  +  Vr,  V»i)  and  equation  (4. 

*  %  “  " 


2) 


takes  the  form 


\-k  --T<7«7f  +  v:-v'>)  +  ^f(y'Z-ri-V ]M4mei 

L  *  **  (4.6) 


where  the  function  f,  expressed  in  terms  of  the  new  variables, 
is  independent  of  the  center  coordinate  R  and 


H 


Equation  (4.6)  and  (4.7)  together  with  the  boundary  conditions 
on  M4  at  z  3  0  determine  the  corresponding  fourth  order  statis¬ 
tics  of  the  field  in  the  random  medium.  To  date,  no  analytical 
solutions  have  been  derived  for  equation  (4.6),  except  for  the 
special  case  H( J  )  =  a  £  which  has  been  solved  by  Furutsu(li) 
using  functional  techniques. 

In  section  2,  we  shall  extend  the  "quasiparticle  distri¬ 
bution  function"  concepts  to  a  4-point  correlation  spectrum, 
i.e.  the  Fourier  transform  of  4-point  correlation  function,  to 
obtain  an  analytical  result  for  the  case  H(  £  )  =  a  £2,  An 
approximate  solution  for  arbitrary  power  spectrum  using  4-point. 
Slow  spectrum  will  be  derived  in  section  3.  In  the  last  section 
of  this  chapter,  we  will  present  a  numerical  scheme  to  solve  the 
fourth  moment  equation  (4.6)  for  an  arbitrary  power  spectrum. 


II.  The  Exact  Solution  for  the  Case  H(  £  )  =  at2 


We  shall  define  the  F  and  S  functions  corresponding  to  the 
fourth  moment  M4,  vizi 
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where  R  and  are  denoted  as  the  "slow  variables",  and  ?  and 
r,  as  "fast  variables".  One  observes  that  F  and  S  form  a  4- 
fold  Fourier  transform  pair  with  respect  to  R,  ra  and  kj,  ,  k,.  . 
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Furthermore,  the  "macroparticle  coordinates"  can  be  obtained 
by  the  following  relations 
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For  the  special  case  ,  H{  £ )  =  a E 2,  one  obtains 
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Substituting  (4.11)  into  (4.6)  and  using  (4.8b),  we  can  derive 
the  defining  equation  for  S  asi 
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(4.12) 


The  characteristic  trajectory  equations  corresponding  to  equa¬ 
tion  (4.12)  are 
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Equation  (4.12)  can  be  described  as  the  slow  quasiparticle 
distribution  function  in  a  two-dimensional  phase  space  spanned 
byk*,*,-,?  ,r,  .  Along  the  trajectories,  the  S  function 

V  w  V  M 

obeys 


S  (*.  ,f  n,  i)  =  S  °)  * 
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Q  °  ~  (4.14) 


which  is  essentially  the  same  expression  as  described  in  the 
2-point  correlation  case  except  that  6-fold  spatial  directions 
are  involved.  For  simplicity,  we  shall  consider  only  the  two 
dimensional  case  which  can  be  extended  to  the  6-dimensional 


case  without  great  difficulty 


For  incident  wave  at  z  =  0  of  the  form 
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the  four  point  correlation  function  can  be  calculated  as 
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Using  (4.8b),  one  obtains 
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where  K*  and  kr|  are  normalized  by  k0  .  Substituting  (4.16) 
into  (4.14)  and  appling  the  trajectory  equations  defined  by 
(4.13),  we  get 
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The  fourth  moment  defined  by  (4.15)  can  be  derived  by  substitu 
ting  (4.17)  into  (4,8b).  For  example,  let  us  consider  a 
special  case  when  rx  =  P  =  r4  *  0,  then 
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and 


C  K,o,o.  o,  l)  = 


O^y 


4-oir, 


0  cj  feft  ci  far, 

(2K)Z 

_  _£i 

6  ^-3 


(4.19) 


4  6 


For  the  ionospheric  propagation  case  described  by  Fig. 1.1, 
the  S  function  obtained  in  equation  (4.17)  will  be  used  to  ob¬ 
tain  the  distribution  function  at  the  bottom  of  the  random  slab 
at  z  =  zx.  The  quasiparticles  will  redistribute  while  propaga¬ 
ting  in  free  space  and  obey  the  following  equation: 


&r,,  r2j  J)**  O 

(4.20) 


The  trajectory  equations  becomes 
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where  zz  is  the  location  of  the  receiver  on  the  ground.  One 
notes  that 
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using  (4.17),  (4.22)  and  (4.21),  one  obtains 
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where  A  z  =  zz  -  zt  is  the  distance  from  the  bottom  of  the  slab 
to  the  ground.  For  the  special  case  ?  =  r,  =  rz  *  0  the  S  func 
tion  can  be  calculated  from  (4.23)  as« 


S  (fta,  &r.  ,  0/  30 


f <-** V-  3/  -  fe* 


(4.24) 
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and 

M4  (  K,  0,  o,  o,  }2) 


where 

A  =  T°lr^+1 X 

E  ,  J_ dr3,*+  -i-  (t-diiS*) +*■(*.**'->' M.26) 
U/  y  o  «r 

tfS,3) 


One  thus  concludes  that  if  the  4-point  statistics  of  the  wave 
field  is  given  at  z  =  0,  the  corresponding  4th  order  statistics 
at  the  ground  can  be  obtained  from  the  z  evolution  of  the  qua¬ 
siparticle  distribution  function  from  equation  (4.14)  and  (4.22). 


III.  Approximate  Solution  for  4th  Moment 

The  defining  equation  for  the  4th  moment  can  be  rewritten 

as  x, 

£  J.  -  £  (  V,  vf  +  Vr,  vri)l-  (H  (rx+  if)t  H  (rx-±f)]j  M+ 

r  -  Jb!  [W  +i?)  1-  H  (r-if)-H  (nmy-Hfr-nj]  (4-27) 

4  L 

m4  . 


-  Ill 


Inserting  (4.8b)  in  (4,27),  one  obtains 


£  JL  +  •  ?r,  *■  -^-[h  l*i  +  if)  t  H  ( 5-  if)]J  S 

.a  f  -iffc.-R  +  Ar.-rr) 

=  I  M4  e  *-  * 


JfZ  d  y; 


(4.28) 


where 


J*  H  (r;  +-£f)  +  H(r,~  if)  -  H  *}J) 

-  H  Cr.-r*) 


since 


H  (?)  =  are  :f<£  («;  fi -  e  * '] 


one  derives  from  (4.28)  after  some  manipulation 

=  4 ^ 4 > [c*<>  4 - >>  - Ctstl.  j 

e *4‘5  d  fe 

and  .  . 

c  -Z  (*.  R  ■<■*<■,  -r.) 

J  e  - -  ~  dR  Jr, 

=  V-s-^  ^JS 

*'-*r,,ra,  j;  J  6' 


(4.29) 


(4.30) 


(4.31) 
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From  (4.28)  and  (4.31)  we  readily  derive  the  integral 
equation  for  S  as» 

{jf  +  T.  •  V1  +  -fe  • v?  +  T-tH  +  H  (y^}s 


J 


i  y  J*' 


32) 


Equation  (4.32)  contains  an  extinction  term  H(rz+  f  )+H(r,  -hf3) 
which  will  reduce  the  quasiparticle  distribution  function~while 
redistributing  it  in  k  R  ,  f  ,k,*  ,  r  space.  The  "scattering 
term"  on  the  right  hand  side  of  the  equation  describes  the  con¬ 
tribution  to  the  total  S  function  at  the  "momentum  coordinates" 
Kr  «  from  a  "momentum  coordinates"  k*  ,  k'  -  kr  .  The 

t  ^  'V'  ^  ^  1 

complicated  form  of  equation  (4.32)  makes  it  impossible  to  ob¬ 
tain  a  solution  in  analytical  form.  We  shall  consider  first 
the  simplest  case  of  a  plane  wave,  for  which  obviously  V*  M4  =0. 
Vf  M4  can  also  eliminated  from  the  equation,  (  ?  is  only"  a 
parameter)  by  setting  it  to  zero  because  in  this  case,?  =  0, 
the  points  f,  ,  lie  in  the  plane  z  at  the  vertices 

of  a  parallelgram  centered  at  the  point  R  and  with  sides  -  fi  '  = 
-  ?1SS  r2  ;  ?,-?»'  -f2  =  r..  Therefore  equation  (4.27) 

V  V  V  ^  V 

takes  the  form 


(4.33) 


where 


X  (r,.n,s)  =  2H(r.)+2H(n)-H(*+n)-H(yi-n) 

(4.34) 


one  notes  from  equation  (4.5)  that 

r>s  t  tR-  <*'>] 


which  suggests  that  r,  and  ra  can  be  treated  as  ''slow"  and 
"fast"  variables  respectively.  In  this  case,  one  encounters 
the  same  degree  of  difficulty  using  either  the  "slow  spectrum 
S"  or  the  "fast  spectrum  F"  to  solve  equation  (4.33).  In 
order  to  maintain  the  consistency,  we  shall  adopt  the  "slow 
spectrum  procedure"  as  described  previously. 


Following  the  same  steps  as  in  (4. 27)-(4.32) ,  one  obtains 
the  defining  equation  for  S  function  for  the  plane  wave  case  asi 


+  -^H  (}.r,)jS  (**>.*.}) 


=  J  #  (*’->  C1  -  ^  *'•  5  J  S  (*-  *r„  4,  Sj) 


J  ft' 


(4.36) 
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On  the  trajectory  defined  by 


*1 


(4.37) 


equation  (4.36)  becomes 


—  H  Cna,) ^  (.-kr.r^ilS) 

A}  2  ~  ~ 


=  J  (  £r,  ,  *i  ^ 3 y  } 

V  .  AT  ^ 


(4.38) 


where 


Jli„r,(i),3)  =  -A.J  §  (*')[l-ccs  *' ’-nr*)] 

5(Kr^',ni3ii)JA’ 

„  *  '  m,  '  ^ 


(4.39) 


The  solution  of  (4.38)  has  the  form 

5(*r,,nt5>/5;  =  5 <•>,£> j  e 

rl  .  i**  f  (nd^I 

+  )  <?  *' 

o  -  — 


r3 

-  (r.<P)d  J 


/  (4.40) 
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Substituting  (4.37)  into  (4.40),  one  obtains 


2  C  3- 


S  (*r„  r„  3;  =  5  (  fer„  ^  l  ,  O)  ^r,  [-  L’jj 

H  <3  -  Is  Li-*))  J 1  +  (-  ( 

H  1 5  -  %  £3-5  0  4j  J  ( fc, ,  5  -  $£  «- A  O  j <3  ’ 


(4.41) 


Equation  (4.41)  together  with  the  z  =  0  boundary  condition, 
S(o,  r4  ,  kr  )  =  (kr,  )»  gives  the  solution  for  the  S  func- 
tion  at  distance  z.  The  first  approximate  solution  can  be 
obtained  by  substituting  the  z  =  0  value  of  the  S  function 
into  equation  (4.39),  whence  one  obtains 

J  2')  -  ^  §  (*«)  * 

[l-cos  £r,-(r2-fyoQ-}f4'4t 

and  from  equation  (4.41) 

-  %^H(n)  J 

S  (Ar,,  *i,  3)  *  S  (  fer,)  6  + 

-V  -V  **' 

Joe  iJ3  ~  * 


[1  -  cos  ft, r,  a-i'»]  cl  3 ' 


(4.43) 
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One  notes  that  the  Fourier  transform  of  equation  (4.43),  which 
gives  the  approximate  solution  of  M4,  is  not  a  symmetric  func¬ 
tion  of  r,  ,  r x,  .  This  violates  the  symmetry  property  of 

A>> 

Therefore,  it  is  logical  to  consider  the  Fourier  transforma¬ 
tion  with  respect  to  r2  ,  i.e.  via  the  fast  spectrum  procedure. 
Let  us  define 


F  (n,-6r.J)=  M. 


;e 


4V  ** 


(4.44) 


One  obtains  a  solution  for  F  similar  to  that  described  in  (4.43), 
except  for  the  interchange  r2<— *  r,  and  kr«-»kr  .  The  approxi- 
mate  solution  of  M4  can  then  be  written  as 

M4  Cn,r2,|;  =  J  [Js  (K j  r,  + 

.0  (4.45) 

J*  F  ( **•  ^Sx  '  5^  0  A  "  *  J  r2  J 

which  maintains  the  symmetry  character  of  M4  with  respect  to 
r,  t  .  (cf.  Tatarskii,  1) 


IV.  The  Boundary  Conditions  and  Steady  State  Solution  for  Plane 
Wave  Case 

We  note  that  when  the  pairs  of  point  ( f,  ,  '  )  and 

v  %■  v 

move  an  infinite  distance  apart.  It  is  obvious  that  the  fields 
corresponding  to  these  pairs  of  points  become  statistically  in¬ 
dependent,  thus  one  has 


*  rt 
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at  I  Yx  1  — ►  oo 

M4  (r,,ra>  l) 

%» 


<  $  <  3.  t, )  a &)f  «,«')> 

M2(S.<j)  M/o.n) 


=  I  Mj  (5.H) 


(4.46) 


similary,  we  have 

at  |  r,  l  — *  oo 


M.  («,  r,  .  3)  =  l  Mi  <3*>S)|2 

(4.47) 


we  now  consider  the  properties  of  the  function  f(r,  ,  r*  ,  z)  in 

^  V 

equation  (4.34).  One  observes  that  f(r,  ,rx  z)  is  positive 
everywhere  except  at  r4  =  0  or  ra  =  0,  f(r,  ,  r,  )  =  0,  We 

**  v  <w 

therefore  consider  the  term  K#  /4  f(rt  ,rt  ,  z)  as  a  decaying 
term  which  tends  to  make  the  function  M4  decay,  while  the  2nd 
term  in  equation  (4.33)  is  a  diffraction  term  which  will  diff¬ 
use  the  value  of  M4  among  the  transverse  coordinates.  With 
the  diffraction  term  only,  equation  (4.33)  can  be  written  as 


jm4 

a  I 


•  9  5 

x  ar»  an 

m 


o 


(4.48) 
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the  solution  of  (4.48)  can  be  written  as 

(r,-n';Qr>-rz  ) 

M4  (*■/,<  o;e 

~  ~  (4.49) 

J  v,'d  rz' 

V  V 

Thus  at  z>0  the  value  of  M4  will  be  redistributed  as  descri¬ 
ed  by  the  above  equation. 

One  also  observes  that  due  to  the  decay  term,  M4  will 

decay  to  zero  when  z  >>  1  for  all  r,  and  r,  except  near 

v* 

r 1  =0,  rx  =  0.  However,  because  of  the  diffraction  term  the 
value  of  M4  near  r,  =  0  or  r2=  0  will  be  redistributed  to  all 
other  points  (r(  ,  r2  )  where  f(r,  ,  r2  )  =  0.  Therefore  one 

•w  v  <v 

observes  that  M4  will  decay  to  zero  for  large  distance  z  unless 
M4  is  of  the  form  (cf.  Lee  and  Jokipii,  68,69.) 

(4.50) 


M4(»;,rJ,3;  = 

A*-  V 


7 


For  this  range  the  diffraction  term  is  zero  and  M4  will  reach  a 
asymptotic  solution.  In  order  to  satisfy  the  boundary  condi¬ 
tions  in  equation  (4.47),  we  must  have 

at  I  Vj  |  — * 

M+  *  |  M2  Wi  (2,ao)  +  'l”z(  2,  + 

(4.51) 
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therefore 


=  in,  (»,*)  +  -*• 

-  7n3  (23  -f  Tn3  (§3 

=  7n,CI>V7^+  l  M2  -Wi 


(4.52) 


at  I  J* I  — *  00 

M4=  |  M2  (5-^)|2=  w,v'3-n)+  Vi2t3,~) 

1 1«3  (3) 


(4.53) 


therefore 


m,an>=  I  M2Cl.r)lZ-yn2^^)-Wi(V 


m,(j.->-  I  M*  (3,~>| -Wa  <*•••>- 


(4.54) 


Substituting  (4.54)  into  (4.52),  one  obtains 

M4  (J»l)=  |Mz(S.r.)|  +|M2a,0)| 

-  |  Mi 


(4.55) 


For  the  plane  wave  case,  we  obtain  the  following  expression 


M.  U»i,  *.r,)  = 

/V  V 


-*£[»<*  h  -¥&(r‘U 

0  +  £  (4.56) 

-ILL/HO 

e 


After  the  wave  leaves  the  random  slab  with  thickness  zA  and 
propagates  in  the  free  space  region,  we  have  the  differential 
equations 


M4  =o 

9  <)  Ko  ~  ~ 


(4.57) 


the  boundary  conditions  for  M4  are  given  by  equations  (4.46). 
However,  M2  appearing  in  equation  (4.46)  now  takes  the  form 


_ 


Mz(i ,r)=  M 2(ll9r)  e 


kojj 

2 


[A(t>) -A 


(4.58) 


since  in  free  space  M2(z,r)  is  unchanged. 

The  solution  of  (4.57)  can  be  written  immediately  in 
terms  of  the  value  of  M^z,  r^  ,  r2  )s 
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m4  ( 5;  =  in'.  5’.  5; 


-x  4s- 


j  n'J  ^ 


,  (4.59) 


for  z>0  ,  the  value  of  M4  is  redistributed  among  the  trans¬ 
verse  coordinates  until  it  reaches  a  steady  state,  if  such 
exists.  One  finds  that  at  steady  state,  M4  obeys 


M4  =  |  M2  ( 5,  ,  |2+  j  Mz  (§i<  >p) 

-  |  M2  <‘,•60, 


The  general  solution  of  M4  is  too  complicated  for  compu¬ 
tation.  In  order  to  simplify  the  problem  and  still  keep  the 
amin  features  of  the  solution,  we  shall  take  the  four  points 
f,  ,  fz ,  f,' ,  on  the  z  =  constant  plane  along  a  straight  line 

such  that  r,  and  ra  in  equation  (4.33)  become  scalar.  Intro- 
ducing  the  dimensionless  variables 


3 


r,  = 


-  y'z 

^  *  —cr 


(4.61) 


where  $,  is  a  parameter  related  to  the  scale  size  of  the  irre¬ 
gularities.  When  written  in  terms  of  these  variables,  eq.(4.33) 
becomes 


JL 

an 


+  -j-  (o)fjMA  =  0 

(4.62) 


where 

f  (n  .n)  =  2  -  ~  [2  A  (r.)  +  2  A  (n)  -  A  (r.r^) 


3  l 

one  observes  that  =  k  kpj^A(O)  is  a  dimensionless  quantity. 
For  convenience,  we  shall  drop  all  the  "bars"  and  rewrite  (4.62) 
as 


54  _:_L.  _1_ 

is} 


vf  (f,,  n  (r,,n,z)  =  0 

J  *  (4.64) 


Since  M4  is  a  complex  field,  it  is  convenient  to  divide  it  into 
real  and  imaginary  components  by  writing  M4  =  Mr  ,  one 

obtains  two  coupled  partially  differential  equations 
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and 


A  («)  »  2jni  -fen'  <:  V > 


(4.68) 


however,  in  situ  measurements  indicate  a  more  realistic  form 
is  given  by  the  modified  power  law  spectrum  defined  through 


(4.69) 


with  Jclo<<  1%.  This  spectrum  is  flat  for  k  <  kL„  »  is  a 

power  law  with  index  -p  for  kL  < k  <  *m  »  and  is  exponent i- 

0 

cally  decay  for  k  >)<„,.  L0  =  IAl,  is  the  outer  scale  and  10= 
1/lCjjj  is  the  inner  scale  of  the  irregularities.  In  the  iono¬ 
sphere  case,  usually  2<  p<4  and  p  =  11/3  corresponds  to  the 
modified  Kolmogorov  spectrum  (Von-Karmann  Spectrum). 


For  the  Kolmogorov  spectrum,  the  power  spectrum  can  be 
written  as 

ft* 

§  «y=  o  o33Ct2L0/*(l+  ^Li)  Vte  *” 

(4.70) 


1-1  <  v> 


with 
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if  L#»l0  one  obtains 


Hi?) 


A  Co)  -  A  (?) 

<£2>  Lo 


H 


K’a  (?) 


3- 


71) 


A  (o)  So-315  <£*>  Lc 

(4.72) 


where  denotes  a  modified  Bessel  function  of  the  second 

kind. 

If  we  choose,  Ld  =  1  =  300m,  and  <ez>  * 
the  value  of  P  =  ti  l2<?>  ranges  from  about  8.5  for  f/fp  *  100 
to  sbout  86.5  for  f/fp  =  10.  The  value  of  •y*  in  eq.  (4.65)  is 
equal  to  0.0987  £  for  the  Von  Karmann  spectrum,  0.16  ^  for  the 
power  law  spectrum  with  p  =  4,  and  0.2821  for  the  Gaussian 
spectrum. 

We  shall  adopt  the  power  spectrum  given  by  equation  (4.69) 
with  p  -  4  and  present  a  numerical  solution  for  the  fourth  mo- 

__  £  f 

ment.  The  value  of  °r  =  !sK0Ls,A(O)  has  been  chosen  to  be  3.5. 
Fig.41  depicts  the  normalized  function  f(rt  ,  r*  ) ,  defined  by 
equation  (4.63),  as  a  function  of  r,  ,  r3  .  One  recalls  that 
f(rt  ,z 2.  )  represents  a  decaying  factor  which  is  positive  every¬ 
where  and  equal  to  zero  when  r,  and/or  rt  =  0.  The  initial 
field  M4(r,  ,  r4  ,  0)  =  1  is  not  shown  in  the  figures.  The  evolu¬ 
tion  of  the  4th  moment  at  z  =  0.1125,  0.1725  and  0.24  is  shown 
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VON  KARNANN  SPECTRUM  F<R1,R2> 

■X 


F<R1.R2> 


Fig.  41.  f (r.  ,  r,  )  vs.  r*  and  r,  5  Power  spectrum  with 
p  =  4. 
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Fig.  42.  M4(rj  ,  »  z)  vs.  rx  and  r*  at  z  =0.1125 

a  plane  wave  propagating  in  a  random  medium  with  a 
power  law  spectrum,  p  =  4,  T  *  3*5. 


in  Fig. 42  -  Fig. 44.  One  observes  that  the  field  decays  every¬ 
where  except  near  the  neighborhood  of  rt  and/or  r2  =  0.  The 
redistribution  phenomena,  caused  by  the  2nd  term  (diffraction 
term)  in  equation  (4.64),  forces  the  points  near  r*  and/or 
r2  =  0  "move”  to  the  decaying  region.  Therefore,  if  the  evo¬ 
lution  continuous  to  a  large  distance,  the  final  field  will 
tend  to  reach  a  steady  state  as  shown  by  equation  (4.55). 
However,  we  are  not  interested  very  large  distance  propaga¬ 
tion  in  this  paper.  The  extend  of  the  normalized  in  random 
slab  will  be  choosen  such  that  L  =  0.24  and  take  the  solution 
will  be  evoluated  at  z  »  z',  whence  the  wave  will  continue  to 
propagate  from  the  bottom  of  the  slab  to  the  ground.  Fig. 45 
48  depict  the  free  space  propagation  at  z  =  0.4275,  0.6525, 
0.8775  and  1.1025.  The  diffraction  effect,  resulting  from 
the  interference  of  the  distorted  wave  front,  tends  to  focus 
and  defocus  the  field.  This  focusing  and  defocusing  phenome¬ 
na  will  only  occur  for  strong  turbulence  case,  which  will  be 
shown  clearer  in  next  chapter. 

VI.  Discussion  and  Conclusion 

The  4  point  slow  quasiparticle  distribution  function, 
defined  as  the  6  dimensional  Fourier  transform  of  the  4  point 
correlation  function  w.r.t,  the  slow  coordinates  R  and  rj,  is 
used  to  obtain  the  4  point  wave  statistice.  Pictures  of  the 
slow  q.p.  flow  are  essentially  the  same  as  those  in  the  2 
point  correlation  case  except  we  are  now  dealing  with  a  four 
dimensional  phase  space  kR  ,  R,  k».  ,  r,  .  An  individual  slow 
q.p.  moves  along  the  trajectory  defined  by  eq.  (4.13)  with 
"momentum"  kR  in  the  f  direction  and  ky,  in  the  r4  direc¬ 
tion.  The  structure  function  is,  in  general,  dependent  not 
noly  on  the  "coordinates"  f  and  r2  but  also  dependent  on  the 
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Fig.  47.  M/j  vs.  r.  and  r z  a* 
bottom  of  the  slab. 
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Fig.  48.  M4  vs.  rt  and  at 
from  the  bottom  of  the  slab. 


"coordinate”  r,  associate  with  the  "momentum"  .  For  the 

special  case,  H(  £  )  =  a  E2,  one  finds  that  the  structure  func¬ 
tion  is  dependent  only  on  f  ,  and  of  the  form  f(rv  ,  r,  ,  f  )  = 

a  f\ 


Knowledge  of  the  S  function  will  provide  all  of  the  4 
point  field  statistics.  For  example,  the  second  moment  of  the 
field  intensity  can  be  obtained  as  follows* 


<  iZ>  =  M4  (o,o,o,o.$)  = 


J 


(4.73) 


Therefore,  the  second  moment  of  the  intensity  can  be  viewed  as 
the  total  number  of  slow  q.p.’s  located  at  R  =  ra  *  r2  =  J  =  0 


The  2  point  correlation  function  of  the  intensity  at  R 
=  0  can  be  obtained  as  follows* 


^  Ic*oK*0>  s  ^4  (  0  'r<  *  °'  $  ) 

r 

=  s  e 

d  d  fer, 

(ilTj1 


(4.74) 


One  notes  that  when  f  =  r  2  *  R  ■  0,  then  rj=  x^-  X£®  xj-  X2 ' 

xi’  -  x2. 


The  Scintillation  Index  S4,  defined  as  the  normalized 
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a ^  t ^  v2  y  \  jt 

variance  of  intensity,  i.e.  S4  =  4.1  >"<J-->Al>,  can  then  be 
obtained  from  the  S  function  using  eqs.  (4.73)  and  (4.74).  As 
an  illustration,  we  shall  calculate  the  1st  and  2nd  moments  of 
intensity  and  the  associated  S4  for  the  case  H(  £  )  =  a  £2,  using 
the  S  function.  Fig.  49  and  50  depict  ^I>,  and  the  corres¬ 
ponding  S4  for  a  collimated  beam  propagating  in  a  random  medium 
characterized  by  H ( £  )  =  a  t1  with  a  =  1  and  a  =  5  respectively. 

A 

(In  the  following  figures,  curve  1  stands  for<l>,  curve  2  for 
<Ia>and  curve  3  for  S4.)  One  observes  that  both  <I>  and  <X*> 
decrease  as  z  increases,  while  the  scintillation  index  increa¬ 
ses  monotonically  with  z.  The  case  of  a  focused  beam  propagat¬ 
ing  in  a  random  medium  is  represented  in  Figs.  51,  52  with  dr~ 
oil  =  1,  a  =  1  and  5  respectively.  Focusing  and  defocusing 

A  *2 

effects  in  <I> and  <1  >are  clearly  revealed  in  these  figures. 

The  scintillation  index  is  greater  in  the  case  of  collimated 
beams  in  the  range  z  <  2,  approximately;  for  z£2,  focused  beams 
display  a  smaller  scintillation  index  than  the  collimated  beams. 
One  thus  concludes  that  for  short  distance  propagation,  a  colli¬ 
mated  beam  is  to  be  per f erred  while  a  focused  beam  is  perfer- 
able  for  long  distance  wave  propagation  in  a  random  medium. 

In  ionospheric  propagation,  we  are  interested  in  finding 
the  received  signal  intensity  and  its  associated  S4  at  ground 
level.  A  deterministic  wave  signal  is  randomly  modified  by  the 
ionosphere  and  when  the  resulting  modulated  wave  front  enters 
the  free  space,  between  a  random  ionospheric  slab  and  the  gr¬ 
ound,  diffraction  modifies  the  wave  statistics.  In  order  to 

A  A  2 

show  this  effect  we  present  results  for  ^I>  ,  4  I  >  and  S4. 
Figs.  53  and  54  depict  the  case  dr  =  1,  a  =5  with  the  thick¬ 
ness  of  the  random  slab  being  L  =  0.25  for  collimated  beams 
and  focused  beams  respectively.  One  observes  that  the  focused 
beam  displays  saturation  at  z  =  0.75,  where  z  measures  the  dis- 
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tance  from  the  bottom  of  the  random  slab  to  the  receiver.  For 
z  smaller  than  approximately  1.8,  the  collimated  beam  yields  a 
greater  S4  than  the  focused  beam.  However,  as  z  1.75  the  fo¬ 
cused  beam  displays  a  smaller  S4.  Figs.  55  and  56  for  a  =  10 
display  the  same  effects.  One  obtains  a  higher  scintillation 
index  and  observes  a  saturation  at  z  =  0.75  for  the  focused 
beam  case.  It  is  of  interest  to  plot  S4  vs.  z  for  different 
turbulent  levels  in  the  ionospheric  slab.  Figs.  57  and  58  de¬ 
pict  34  vs  z  for  a  =  2  to  16,  in  steps  of  2,  for  the  collimated 
beam  and  the  focused  beam,  respectively.  One  observes  that, 
for  the  focused  beam,  S4  always  saturates  at  z  =  0.75.  At  z  = 
1.75  the  focused  beam  start  to  display  a  higher  scintillation 
index  than  the  collimated  beam. 

When  the  slab  thickness  increases  to  0.5,  the  random  wave 
fields  emerging  from  the  bottom  of  the  ionosphere  are  increas¬ 
ingly  more  turbulent.  The  corresponding  S4  in  this  case  is  high¬ 
er  than  for  the  case  L  =  0.25.  Figs.  59-62  depict  ^I>  »^IJ.> 

and  S4  vs.  z  for  a  =  5  and  10  respectively  ;  one  observes  that 
the  focused  beam  saturates  at  z  =  0.6  approximately.  The  focu¬ 
sed  beam  has  a  smaller  S4  than  the  collimated  beam  in  the  range 
z  greater  than  1.5.  The  influence  of  the  turbulent  strength  a 
on  the  scintillation  index  is  shown  in  Fig.  63  and  Fig.  64  for 
a  collimated  beam  and  a  focused  beam  ,  respectively.  The  value 
of  a  ranges  from  2  to  16  in  steps  of  2j  one  notes  that  at  z  = 

0.6  the  scintillation  index  of  a  focused  beam  saturates  for  all 
values  of  a.  At  a  distance  z  =  1.5  the  focused  beam  starts 
to  have  a  smaller  S4  for  all  values  of  a. 

The  wave  front  of  the  wave  emerging  from  the  bottom  of 
the  random  slab  is  greatly  dependent  on  slab  thickness.  As  a 
comparison,  we  present  calculations  of  S4  vs.  z  for  different 
values  of  thickness  ranging  from  L  =  0,125  to  1  with  increments 
of  0,125  as  depicted  in  Fig.  65  and  Fig.  66.  One  observes  that 
for  the  focused  beam  case  the  saturation  distance  decreases  as 
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Fig,  56,  a  modulated  focused  beam  propagating 
in  free  space,  L  =  0.25,  a  =  10,  =  o/.-  =  i. 
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L  increases.  The  distance  at  which  the  focused  beam  displays 
a  greater  S4  than  the  collimated  beam  decreases  as  L  increases. 

The  approximate  solution  for  the  4th  moment  M4  is  pre¬ 
sented  in  Section  III.  In  (4.43),  we  have  replaced  the  true 
slow  q.p.  distribution  function  S(  z,  r  ^ ,  kj.-  k' )  on  the 
right  in  (4.41)  by  the  initial  unattenuated  q.p.  distribution 
function  S(  0,  r^  ,  kr-  k').  It  is  expected  that  the  solution 
so  obtained  will  give  an  exaggerated  result  for  the  intensity 
fluctuation.  To  illustrate  this  fact,  we  shall  calculate  the 
scintillation  index  for  the  plane  wave  case  and  compare  it 
with  the  exact  solution  obtained  numerically.  The  random  slab 
will  be  characterized  by  a  power  law  spectrum  with  p  =  4.  Fig. 
67  depicts  the  scintillation  index  calculated  from  the  approxi¬ 
mate  solution  (4.43)  for  the  cases  L  =  0.06,  0.12,  0.24.  The 
general  behavior  of  scintillation  index  agrees  with  its  well 
known  feature.  Fig.  68  shows  the  comparison  of  the  approxima¬ 
te  solution  (solid  line)  with  the  exact  solution  (dash  line). 

As  expected,  one  observes  an  overestimated  intensity  fluctua¬ 
tion  in  the  approximate  solution. 

Boundary  conditions  and  the  steady  state  solution  are 
discussed  in  Section  4.  We  have  applied  these  boundary  con¬ 
ditions  in  the  numerical  evaluation  of  the  4th  moment  as  pre¬ 
sented  in  Section  5.  The  general  solution  of  S4  is  too  compli¬ 
cated  for  computation.  In  section  5,  we  computed  the  solution 
of  M4  for  a  special  case  wherein  the  4  points  x,  ,  x4’,  xz  ,  xx’ 
lie  along  a  straight  line  on  the  z  =  constant  plane.  The  nu¬ 
merical  evaluation  of  M4,  for  a  plane  wave  propagating  thru 
a  random  slab  with  L  *  0.24  and  characterized  by  a  power  law 
spectrum  with  p  =  4,  is  presented  in  Section  5  for  z  =  0.4275, 
0.6525,  0.8775  and  1.1025(Fig.  45-48).  One  observes  that  as 
the  field  evolution  is  continuous  by  extrapolated  to  a  large 
distance,  the  final  field  will  tend  to  reach  a  steady  state  as 
shown  by  eq.  (4.55), 


Fig.  67.  S4  vs.  2  calculated  from  the  approx, 
solution  (4.43 )j  (1)  L  =  0.06;  (2)  L  =  0.12; 
(3)  L  =  0.24. 


Fig  68.  S4  vs.  z,  L  ■  0.24;  solid  line;  approx 
solution;  dotted  linei  exact  solution. 
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Chapter  5  THE  INTENSITY  FLUCTUATION  AND  SCINTILLATION 
INDEX 


I.  Definition  of  Intensity  Correlation  Function  and  Scinti¬ 
llation  Index 


If  the  new  variables  (4-5)  are  introduces,  the  function 
M4 (r^  ,  rz ,  z)  is  expressed  in  terms  of  the  wave  fields  by 
the  formula  (for  f  =  0  ) 


M4  ( Vi ,  ,  l ) 


<$ 


*  * 

if  ( R  + 


Z  ' 


J)$(*-5 T?.i) 


The  2nd  moment  intensity  function  is  defined  by 


~  ~  (5.2) 


which  is  identical  to  the  fourth  moment  M4  when 
or  rz  *  0. 


One  also  recalls  that 


146 


V  »  V  V  V  ‘  v  / 


(5.3) 


M zLl,*,o)  =  <  Id,* )> 


(5.4) 


where 
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correlation  function  as 


One  defines  an  intensity 


Br  (l ?•  ,fi)  =  < i(3.r.)i(3.t,)> -<ia.vxz(3,f,j> 


(5.5) 


and  one  observes  from  equations  (4.5)  that  r  z  =0,  ftsR+JL 

fx*  R-*5-  »  and  x  =  R.  Therefore  ~  ~  2  ’ 

»  *  *  » 
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For  the  plane  wave  case,  one  takes  into  account  that  M2  is 
independent  of  the  slow  coordinate  and  write 


Bj(3,f.,?a)=  M4(S,r,,0)  -  M/  (3,o,o) 

V  v  ^ 


=  M4(5,r.  ,0)  -  <I>2 


(5.7) 


Thus,  M4(z,  r^ ,  0)  enables  one  to  determine  the  correlation 
function  Bi(z,  f,  ,  f,. )  =B(z,  r.). 

^  V  /V 

The  scintillation  index  S4  is  defined  as 


A  .  A  _I 


c  x  <  (i-<f»‘>  M4(i,o,0;_  <f>z 

^4  /  ^  v  2  ”  a — : - 


<i> 


<  i>* 


(5.8) 


which  plays  an  important  role  in  the  scintillation  problem 
and  is  used  as  a  measure  of  fluctuation  level.  For  the  case 
H(£)  a  a  J1,  one  obtains  from  (5.8)  and  (4.25) 

_  _  _r! 

m4(r, 0,0,0, =  ■j^rfra  e  46 


--==e  4P 

2 


(5.9) 
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where 

g  =  (l-oiA-^2)  +  a*($J2 -izb,  +  4>0 

(5.10) 

p=  +  (1  -<**30  +  ^  (2i3z-2*3,V4^) 


the  scintillation  index  can  be  calculate  from  (5.8)  as 

,  A  ,  A  7 

<1  >  -  <I> 
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<■  x  >2 

M2z  (0,  0,  30 


(5.11) 
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II.  Numerical  Solution  of  Scintillation  Index  for  Arbitrary 
Power  Spectrum 

In  our  computation,  we  shall  first  choose  a  Gaussian 
power  spectrum.  The  corresponding  correlation  function  is 
also  Gaussian.  The  main  reason  for  choosing  such  a  spectrum 
is  to  check  with  the  results  obtained  from  the  thin  phase 
screen  approximation  (4  ).  For  the  Gaussian  spectrum,  the 
f(r^  ,r ^)  in  equation  (4.63)  is  given  by 
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rO  =  2  -  2  [ 


(5.12) 


Using  (5.12)  in  the  numerical  calculation  to  obtain  the  scin¬ 
tillation  index,  we  shall  present  some  computational  results. 
Fig. 69  depicts  the  behavior  of  the  scintillation  index  as  a 
function  of  z  (from  the  bottom  of  slab  to  the  receiver)  for 
the  case  =  10.  Curves  1,2,3  correspond  to  slab  thickness 
L  =  0.0<*  0.  12,  0.24  respectively.  For  L  -  0.06,  corresponds 
to  a  thin  phase  screen,  we  observe  that  S4  increases  monoto- 
nically  with  z,  and  finally  reaches  a  constant  value.  The 
results  agree  very  well  with  the  thin  phase  screen  theory. 

As  L  increases,  corresponding  to  a  thick  random  slab,  the 
scintillation  index  grows.  One  notes  that  a  focusing  pheno¬ 
menon  occurs  when  L  =  0.24,  the  scintillation  index  reaches 
a  maximum  before  it  settles  to  some  saturation  level. 

It  has  been  indicated  that,  in  the  thin  phase  screen  th¬ 
eory,  focusing  occurs  only  in  the  strong  turbulence  case(rms 
phase  fluctuation  in  our  case ) .  Qur  results  shows  the 

same  behavior  and  should  be  better,  since  we  take  into  ac¬ 
count  the  thickness  of  the  slab. 

Scintillation  index  curves  as  a  function  of  z  for  the 
case  L  =  0.06  with  different  values  of  *T  ,  equal  to  5,  10, 

30  are  shown  in  Fig. 70.  The  curve  on  the  top  indicates  that 
a  focusing  behavior  appears  for  strong  turbulence  case  even 
though  the  random  slab  is  very  thin. 


The  following  types  of  power  spectrum  will  also  be  con 
sidered 

(a)  Power  law  spectrum  with  power  index  p  =  4 
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below  the  slabe  and  become  smaller  at  large  distance.  Fig.  76 
shows  the  scintillation  index  vs.  for  case  (b)  with  L  =  0.12 
and  y  =  5  (lower  curve),  *JT  =  15  (higher  curve).  The  distinc¬ 
tion  between  strong  and  weak  turbulence  is  obviously  related 
by  the  focusing  phenomenon.  The  Von-Karmann  spectrum  with  an 
addition  Gaussian  bump  greatly  changes  the  behavior  of  scin¬ 
tillation  index  and  is  shown  in  Fig. 77.  Curve  1  depicts  the 
case  k  =  3 ,T-  5,  L  =  0.12  and  curve  2  displays  the  case  k  =  2, 
y  =  5,  L  =  0.12.  One  concludes  that  the  location  of  the  Gau¬ 
ssian  biamp  has  a  significant  effect  on  the  scintillation  index. 
The  focusing  behavior  occurs  more  than  once  in  both  curves, 
which  indicates  the  additive  bump  increases  the  distortion  of 
the  wavefront  emerging  from  the  bottom  of  the  slab.  Fig. 78 
cc^ipares  the  different  behavior  of  the  scintillation  index 
for  the  case  (b)  and  case  (c)  and  for  the  parameters  y  =  15, 
k  =  2,  L  =  0.12.  Focusing  occurs  in  both  curves.  At  , 

spectrum  (c)  has  a  larger  S4  than  spectrum  (b).  Fig. 79  de¬ 
picts  the  behavior  of  S4  for  case  (c)  with  k  =  2  and  TT=  5,  15 
respectively.  One  notes  that  even  for  5,  the  scintilla¬ 

tion  index  vs.  z  shows  a  focusing  phenomenon  which  is  not 

observed  in  case  (b)  with  *)T=  5.  This  apparently  indicates 

.  '* 

that  the  additive  Gaussian  bump  causes  a  strong  scattering 

effect.  In  order  to  show  the  effect  of  the  power  index  P»  we 
observe,  for  the  case  L  =  0,24,  k/  35.44,  the  scintilla¬ 

tion  index  vs.  z  for  p  =  4  and  p  =  11/3  as  displayed  in  Fig.  80, 
the  solid  line  being  for  p  =  4  and  the  dash  line  for  p  =  11/3. 
One  finds  that  the  scintillation  index  vs.  z  has  the  same  cha¬ 
racter  in  both  cases  but  increases  with  the  larger  power  index. 
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Fig.  79.  S4  vs.  z  for  Von  Karmann  spectrum  with 
a  Gaussian  bump,  £  =  2 ,  OT  =  5,  L  *  0.12  ,  (curve  1) 
£  =  2,*T~  15,  L  =  0.12,  (curve  2). 
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Fig.  80.  S4  vs.  z  for  power  law  spectrum; 

(1)  L  =  0.24,  =  8.86,  p  =  4;  (2)  L  =  0.24, 

*  8.86,  p  3  11/3 . 


III.  Application  and  Discussions 

A  power  spectrum,  obtained  from  the  Atmospheric  Explorer 
-  E  (AE-E)  satellite  in  situ  measurements  (70),  is  used  to 
calculate  the  corresponding  S4  for  the  night  time  equatorial 
scintillation.  The  variance  of  N/N  is  estimated  to  be  0.207. 

We  have  normalized  the  power  spectrum  such  that  the  integral 
of  the  fluctuation  power_J^(k)  dk  over  the  observed  frequency 
range  is  equal  to  /(N/N)“>.  For  computational  purposes,  we 
assume  a  flat  spectrum  for  irregularity  wavelength  greater 
than  5  km  and  truncate  the  high  frequency  tail  for  irregularity 
scale's  smaller  than  0.4  km.  The  resulting  power  spectrum  is 
depicted  in  Fig.  81.  In  order  to  compare  the  theoretical  model 
with  the  experimental  results,  we  have  chosen  an  incident  wave 
frequency  137  MHz,  f„0  =  .9  MHz  and  1  =  600  m  such  that  the 
parameter  ajf*=  1/4  k^  12A(0)  in  equation  (4.65)  is  equal  to 
52.73.  The  propagation  distance  z  is  normalized  to  *io3  km  ; 
i.e.  measures  distance  in  102  km  . 

Radar  observations  suggest  that  the  ionosphere  can  be 
replaced  by  a  random  slab  with  an  effective  thickness  -'200  km. 
The  distance  from  the  bottom  of  the  random  slab  to  the  ground 
is  taken  to  be  about  300  km.  Applying  these  parameters  to  the 
defining  equation  for  the  fourth  moment  and  using  the  numerical 
scheme  described  in  Appendix  6,  one  obtains  a  plot  of  the  scin¬ 
tillation  index  vs.  z  as  shown  in  Fig.  82.  A  scintillation 
index  as  high  as  0.35  is  predicted  at  the  bottom  of  the  random 
slab.  As  the  wave  leaves  the  random  slab  and  propagates  in 
free  space,  diffraction  phenomena  develop.  The  scintillation 
index  increases,  passing  through  a  maximum  value  at  a  distance 
<'»250  km  from  the  bottom  of  the  slab,  and  then  saturates.  At 
ground  (  z  =  300  km  from  the  bottom  of  the  slab),  one  obtains 
a  scintillation  index  of  1.14,  which  corresponds  to  a  peak  to 
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peak  fluctuations  of  approximately  32  dB. 

Experimental  results  obtained  by  S.  Basu  et  al.  indicate 
a  scintillation  index  of  0.6  at  137  MHz,  which  is  lower  than 
our  analytical  results.  For  a  density  fluctuation  of  20%  in 
the  equatorial  region,  the  experimental  value,  S4  =  0.6,  seems 
a  little  lower  than  expected.  In  private  communication  from 
Dr.  Basu,  he  has  informed  us  that  his  experimental  result  may 
be  too  low  because  of  instrument  saturation. 

Possible  explanations  of  this  discrepancy  might  bet 

(1)  The  dynamic  range  of  VHF  receivers  are  limited  to  about 
16  dB  -  25  dB.  A  scintillation  index  as  high  as  1.14, 
corresponding  to  peak  to  peak  fluctuations  of  32  dB,  is 
beyond  the  capacity  of  such  receivers. 

(2)  The  in  situ  power  spectrum  was  obtained  at  an  altitude  of 
250  km.  By  extrapolation  to  magnetic  field  conditions  at 
the  scintillation  measurement  site,  Basu  et  al.  suggest 
an  ionospheric  structure  ^ 200  km  in  extent  about  300  km 
above  the  earth.  The  accuracy  of  this  model  assumption 
is  not  evident. 

(3)  In  our  theoretical  work,  we  have  assumed  a  flat  spectrum 

for  5  km  and  have  truncated  the  tail  of  the  power 

spectrum  for  X|<  0.4  km.  This  assumed  power  spectrum  is 
also  a  possible  source  for  -he  above  discrepancy. 

(4)  From  the  investigation  in  section  II,  one  concludes  that 
S4  increases  as  the  thickness  of  the  random  slab  increases 
In  the  numerical  evaluation,  we  have  assumed  the  thickness 
of  the  random  slab  to  be  ~200  km  which  may  be  larger  than 
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the  real  effective  thickness  of  the  ionospheric  slab. 

In  spite  of  the  above  discrepancy,  one  notes  that  the 
numerical  results  for  the  case  of  a  Gaussian  power  spectrum 
are  in  good  agreement  with  the  results  obtained  by  Mercier 
(4).  For  the  case  that  L  =  0.06  which  corresponds  to  the 
thin  phase  screen  case,  one  obtains  a  saturated  S4  which 
agrees  with  the  value,  S4  =  1  -  expi-Z^f),  obtained  by  Mer¬ 
cier  (4).  Our  numerical  scheme  provides  better  agreement 
with  thin  phase  screen  theory  than  the  scheme  proposed  by 
Liu  et  al.  (60).  Therefore,  the  proposed  propagation  model 
and  numerical  scheme  should  also  provide  satisfactory  results 
for  the  case  of  an  arbitrary  power  spectrum.  The  only  real 
difficulties  are  proper  location  of  the  random  slab  and  es¬ 
timation  of  its  effective  thickness. 

For  further  study,  we  shall  apply  the  theory  and  numeri¬ 
cal  scheme  presented  herein  to  the  multiple  frequency,  multi¬ 
dimensional,  inhomogeneous,  etc.  cases.  Consideration  of  an 
anisotropic  power  spectrum  and  of  regime  represent  nonlinear 
effects  in  the  strong  scattering  areas  for  future  investiga¬ 
tion  and  may  provide  further  insights  for  the  theory  of  iono¬ 
spheric  scintillations. 


APPENDIX  1 


Solution  of  (  La-  V  )  f  =  0  where  L 0  is  the  an  unperturbed 
operator  and  V  is  a  random  perturbation  operator  with  <  V  >  =0. 

To  solve  this  equation,  it  is  convenient  to  define  a  stochastic 
Green’s  function* 

(L„-  v  )  Q 


=  1 


(A3  .1) 


For  example,  in  r,  t  space 

r 

A  A 

**  - ^ 


$(r,  t,  r1,  t‘)  f  ( r ,  t‘)  Jr'dt' 


-»  b  c r- S  u-i1) 


A 

For  uniqueness,  the  domain  of  G  will  be  defined  by  the  require¬ 
ment* 


$  .r',  t')  =  0  f.r  t<  t' 


It  is  convenient  to  introduce  the  coherent  unperturbed  operator 
G  0  ,  defined  by 


<  o 
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The  perturbation  operator  in  equation  (Al.l)  can  then  be  treat¬ 
ed  as  the  scattering  part  of  the  propagation. 

By  equation  (Al.l),  we  have  adjoint  operators,  Le  ,  V  , 

^  ,  which  are  defined  by 


(L.*-v*)  <T  =  6  (  Lo-  v) 


(A1.3) 


Therefore, 


[L0-v)  Gr  =  G  (U-V)  =  1 


(A1.4) 


Similarly,  we  have 


Lo  G'o  ~  L0  -  1 

~  (A1.5) 


Multiplying  (A1.4)  by  G  and  using  (A1.5),  we  obtain 

<J.  (  U  -  V  )  G  =  <5-. 


S  ince  G  e  L  o  s  1 ,  we  have 


i 
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Q  =  Q0  +  V  Gr  9  &0  T 


(ai.6) 


where  we  have  defined 


-  £ 

v  G  = 


T  Grc 


(A1.6a) 


A 

The  operator  T  represents  multiple  scattering  by  the  stochas¬ 
tic  perturbation  V  in  the  background  G0  .  By  using  (A1.6), 
we  obtains 


£  (  I  -  Gr„  V  J  -  Go 


(A1.7) 


i. 


i 


Substituting  (A1.7),  into  (A1.6a),  we  obtain 
A  *  -  — I  ^  ^  ~  \  l 

T  =  V  Cr  Q,  =  V  (1-  GCV  ) 

(A1.8) 


or 

A  <v  <v  V  ^  /'T 

T  s  +  =y+v^v 

(A1.9) 
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Note  that  the  inverse  operator  can  be  expanded  in  a  Neumann 
series: 


(1  1  +  5.  V  +  <k.  V  <k.  V  t . 

(A1.10) 


but  it  displays  secular  divergence  that  limits  its  range  of 
applicability. 

By  using  (A1.6),  (A1.8)  and  (A1.9),  we  obtain  ensemble 
average  relations  as  follows: 


$•  -  Go  +  Qe  Vc  Q  =  Gr0  +  G0~T  G\>  (ai.id 

where 

C,  =  <Q> 

A 

T  »  K  T  >  (ai.12) 

Vc4  =  TO-. 

where  Vc  is  a  smoothed  scattering  operator  in  a  background  G. 
On  using  (Al.ll),  (A1.12),  one  obtains: 
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T  =  VcCrG:'  =  vt  Cl -G.Vc)\.Ga~' 

-I  (A1.13) 

=  Vc  Cl-  G.l 'e) 


or 

T  =  vc  (1+  o-.G. 

=  Vc  (1  +  GrVc) 

=  Vc  +  Vc  G,  Vc 


(A1.14) 


Multiply  (Al.ll)  by  L 0  and  make  use  of  (A1 . 5 ) ,  we  obtain 


L0  Gr  -  Vc  C? 


CLo  -  Vc  )  Gr 

<k  CLo-Vc) 
1 


(A1.15) 


Equation  (A1.15)  provides  a  nonlinear  defining  equation  for  G, 
the  nonlinearity  arising  from  the  nonlinear  dependence  of  Vc 
on  G.  Using  (A1.6)  and  (Al.ll),  we  obtain 


Q  =  Gr  +  (k0f 
T  =  T  -  T 


(A1.16) 


L66 


From  (A1.9),  we  obtain 


.  ^  ,  -V  A  <V 

T  -  V  +  VO  V  - 

(A1.17) 


Using  (Al.il)  and  (A1.12),  we  have 

0ro  “  Gr  C 1  +  Vc  G-) 

-  ( 1  +  <3r  Vc  )  £r 


(A1.18) 


Therefore,  we  can  express  (A1.16)  as  follows 


Q-  =  Q  •+  £r  C  1  +  Vc  (k)  T  (1 1  GrVc)  O’ 


(A1.19) 


£r  =0  +  Gr  Tc  Gr 


(A1.20) 


where  Tc,  the  multiple  scattering  in  the  smooth  background  G, 
is  defined  by 


Tc  =  (l  +  v'cG-f'f  (i  +  G'4)"/ 


(A1.21) 
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On  averaging  (A1.4),  we  obtain 

<(U-v;£  >  ••=  i 

=  <  4  (u-v)> 

Comparing  with  (A1.15),  one  infers 

VcQ-  =  <  VQ  > 

From  (A1 . 20 ) ,  we  have 

<v  ^  _  — / 

Vc  =  <  v  G-  >  G 

sr  <vGr  +  V(kTc(k>Gr 

=  <V(^fc> 


(A1.22) 


(A1.23) 


(A1.24) 


We  would  like  to  observe  the  physical  meaning  O:  (A1.20)  as 

A 

follows  from  (A1.4),  we  can  write  G  in  the  renormalization 
from  as 


(  L.  -  Vc  -  V,  )  Q  •  { 

(A1.25) 


Therefore,  the  scattering  is  now  caused  by  a  stochastic  per- 
turbation  *  V  -  Vc  of  the  coherent  renormalized  back- 
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ground  operator  (  L0  -  Vc  ).  Multiplying  (A1.25)  by  G,  one 
obtains 


A  A 


Gr  =  Cf  +  Q-  Vi  Gj  —  Q  i*  Gr  Tc  Ck 


where  we  have  defined 


A  *  ^  _ 

V,  Cr  =  Te  Q 


On  using  (A1.26),  one  obtains 
~  a  *  r  -/ 

Tc  =  V,  G  & 

=  ViC  G  +  QVt  G  )  Gr 

A  A  *  a 

-  V,  +  V,  &  V, 


/V-  A  *  .  ~  / 

Tc  =  V,  11-  CrV,) 


From  (A1.24),  one  finds  that 


=  <  v,  G  Tc  > 


By  using  (A1.27),  one  obtains 


(A1.26) 


(A1.27) 


(A1.28) 


(A1.29) 


■  a  r 
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A  A  A 


=  v,  t  v(  5  V, 

A  A  A 

=  v,  +  v,  Q  V,  * 


Vc  -  <v,  d  v,  >  +  ■■ 


(A1.30) 


(A1.31) 


Since 


Therefore 


S  V-Vi 


=  V  +  V  (V  V  -  <  VG 

_\7(^<V<qr\/>+ . 


(A1.32) 


V  ^  ^ 


\/c  =<'70*V>+<v'Ct  V>Gr  <v  Qv  > 


(A1.33) 


On  expanding  (A1.21)  and  making  use  of  (A1.20)  in  (A1.17),  we 
obtain 


Tc  ~  (  1  -  Vc  Cf  t  Vc  (q-  Vt  Cr  +  *  * "  )  T 
Cl-  Gr  +  QVc(kVc  - . ^ 


(A1.34) 


17  0  - 


A/ 

T 


V  +  VGV  -  <VGV>  +  V<rTc$V 


-  <V<r?cGrV> 


(A1.35) 


Substituting  (A1.24)  in  (A1.34)  and  making  use  of  (A1.35),  we 
obtain 

fc  =  V  t  V  (jv  -  <  V  G-V  >  +  V  GfV  (k  V  + 

(A1.36) 

\/c  =  <vr$-vr>  +  <  v  GrV  G-V  >  ■+*--• 

(A1.37) 


i 


> 


I 
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APPENDIX  2 

In  order  to  find  a  rapidly  convergent  series  for  Vc  G, 
one  considers  the  following  possible  expansions; 

(1)  Expansion  ai 

VcCr  =  <  vGV>Q  +  <vG'V<kV>$  + 

<  v<*V<i  VbV>Gf  -<v<k<v<iV><*V>Gr  (A2.i> 
—  <vGV?>Gjr<v'6rV'>Gf  . . 

as  described  in  Appendix  1. 

(2)  Expansion  bt 
since 


-  <  v  Q  Tc  > 

(A2.2) 

and 

Tc  =  (1  +  Ve  Crf'f  Cl  +  SK)  ' 

(A2.3) 

one  obtains 

Vc=  (l+ViC,),<VGT>  (i  +  <*Vcf' 


where 
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T  =  V  +  VGV  -<vGV>  +vGTt cGV 
-  <vCrTc  G  V  > 


(A2.5) 


One  multiplies  eq.  (A2.5)  by  VG  and  takes  the  ensemble  average 
of  this  equation.  Using  equation  (A2.3),  one  then  obtains 
(assume  <  y  Jn+1  >  =  O  ) 


<v<t  T  >  =  <vCiV>  +  (1+  vcG)tv6vGT  <5-  V>  (a2 . 6 ) 
(1+ 

From  equations  (A2.5)  and  (A2.3)*  one  observes  that 

<  vG  vGrT  G-  V>  *  <vG  V  G  V  G  V>  t  (1  +  VcG) 

<:  vG  VG  v  G  T  G  vG V>  (ltG  14) (A2  * 7  5 


Therefore,  one  finally  obtains  the  expansion 

Vc  ~  ( 1  +  vcGr)}<vGv  >  (1  +  GVc)  '  + 

(‘ltVcG)Z<vGvG,vGV>(ltGVc)  ~f  (A2 . 8 ) 
(1  t\ZcG)  *<vG  v GvGv Gv G  v >  (i  -+-G \ /c)  *+ 

(3)  Expansion  C« 

From  equation  (A2.3)  and  (A2.5),  one  obtains 

-£*  (l+vcG){v  +  vGv-<vGv>  t  vGTiGV 
-<vGTcGv>)  +  {A2,9) 


<  ; 


f 
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Therefore,  from  equation  (A2.2)  and  (A2.9),  one  hast 

Vc  =  <vGTc>  -  Cl+VcGf* \<v(xV>  + 

<VQ  V<$  V>  -  <vCr<  V  (k  V»  +  (A2.10) 

^VCtVQtTcGt  V J>  —  < V&  (1+  GrVc)  ^ 

One  rewrites  equation  (A2.10)  as 

Vc  *  ^  V  £r  fc  >  =  Cl  +  VfcGO  l{<vG&?1‘<VGvGr&Z}‘ 

'  Cl  +GrVc)  V  (it  VcG)  (a2.id 

(1+  GrV'c)H<Vr^Tc>W<  V<?fi> 

=  (1+  VcGr  )H{<V<kV>t<V&V6rVr>}Gt<*Vc) 

where  C  1  “*  Po  J 

po  *  (  1  +  <  V<T  V<qrTc  £  (!*"  (A2.12) 

<  V<qrTc>~' 

The  term  <VGVG*£GV>  in  eq.  (A2.12)  can  be  written  asi 

<  v<kv<i-Tc$V>  »  (It  Vc&)-'  f<v<rV'4V<Jv> 

t  <v$v'&  v&v<kV  >  -<v^v$<v6j>$VJjw.  i3) 

<l  +  fr\fc)"'  Cl-  Pi]'1 

where 

Pi  »  (1  +  Vc  SO'1  {  <  V  & v<k  V<T  Tc  G?  v  G  V> 

-<v<k$r<k  <  vGTc.  <k  v>  Gr  ?>}  (A2.14) 

(1+  GVC)~'<  vG  vGrfe.<kv>~' 


1 


One  also  assumes  that  n  is  a  pure  stochastic  variable  indep. 
of  r  and  z.  The  defining  equation  for  the  average  Green's 
function  is 


* 


<n( 5-  >  =  1 


(A2 . 16a ) 


or 


•  2  &  l.  v/  r 

x-n  c * 


=  i 


(A2. 16b) 


Taking  the  Fourier  transform  of  (A2.16b),  one  obtains  the 
algebraic  equation i 


(b  +  i*  -  Vc)  f 


=  1 


(A2.17) 


or 


9- 


b  +  ^  6  -  Vc 


(A2.18) 


Mfcwmiiaiiaatttyflir  ti  i  ir 
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In  order  to  evaluate  g,  one  requires  information  about  Vc  . 

The  most  appropriate  expansion  for  V£  g  will  be  determined  via 
numerical  investigation. 

(1)  For  expansion  a,  we  define  the  following  approximations: 


1st  approx. « 
2nd  approx. : 
3rd  approx. : 
4th  approx. : 


Vc$  =  ?<v‘> 

vc|=  r«r>  -*<*'>'} 

Vc ^  =  [2nJ  appro*.}  +  §*  £  O  } 

Vc$  =  {  3rd  Appro*.}  1-  CD  J* 


For  the  case  that  ns-v  is  uniformly  distributed  between  -a 
and  a,  one  finds  that  Vc  g  tends  to  converge  as  shown  by  the 
solid  line  of  Fig.  1.  Note  that: 


1st  approximation 
2nd  approximation 
3rd  approximation 


is  valid  when 

is  valid  when 

is  valid  when 


f<v'>  =  <0.2 

$  o .  ^ 

92<  v\>  <  o-v 


(2)  We  define  the  following  approximations  for  expansion  b: 


1st  approx. 
2nd  approx. 
3rd  approx. 
4th  approx. 


(1  +  Vc$  V*> 

Vc$  =  J  1st  }  -t  CitV'cJj 

Vc %  *  App»i>y.}  +  (1+Vc$) 

VcJ  *  [3^  *PPr6*}  +  (ltVc$j 


e.0  S.00O-01  1.00  iso  a.  ee 

Fig.  A2.2.  Vcg  vs.  expansion  ci 

(1)  1st  approx.}  (2)  2nd  approx.;  (3)  3rd  approx.; 

(4)  4th  approx.;  — 
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If  n  obeys  a  Gaussian  statistics  with  a  distribution  func¬ 
tion  of  the  formi  „  . , 

f  n  i 

A e~  J 

one  finds  that  all  expansions  described  above  lead  to  poor 
convergence.  As  an  example,  we  calculate  the  1st,  2nd  and 
3rd  approximations  for  expansion  c  and  depict  it  in  (Fig.3). 
Let  us  observe  the  exact  solution  of  G  for  equation  (A2.15), 
when  n  is  Gaussian,  we  have 

_  <.»•*>  C3-3' 

1 

As  a  power  series  expansion,  one  writes 


rer  series  expans i 

. J 

If  the  series  is  truncated  after  a  finite  number  of  terms, 
we  have 


AS 


3-3 


The  "secular"  behavior  of  G  in  this  case  explains  somewhat 
why  "moment'  expansions  for  Gaussian  statistics  are  not 
suitable. 


Consider  a  Gaussion  statistics  truncated  at  £.|a|  namely* 
•f(n)s  /\  0  ~  (“E-)  -ck<Tl  < 

The  results  of  numerical  calculations  of  Vc  are  shown  in 
Fig.  4  -  Fig.  6.  One  observes  that  the  convergence  of  the 
series  expansion  is  dependent  on  the  ratio  *'/<**> 
where  <  n*>  denotes  the  variance  of  n.  For  Fig.  4  to  Fig.  6, 
we  use  59796,  3.62336,  3.05054  respectively,  and 

find  that  a  smaller  ratio  gives  a  better  convergence. 


[ 


#.74«-ei 


2  83*-97  L_ _ J _ I _ I _ I 

0.9  7  S0«— 01  1  50  2.25  3.00 

g*  <*> 

Fig.  A2.3,  Vcg  vs.  ga<v>;  n  is  Gaussian. 
Expansion  ci  (1)  1st  approx.} (2)  2nd  approx.; 
/  (3)  3rd  approx.; 


Fig.  A2.4  Vcg  vs.  g**vV;  n  obeys  a  Gaussian  statistics 
truncated  at  ♦  | a |  ;  a*  /  4.  4.598.  (expansion  c) 


***** 
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APPENDIX  3 


The  limit  of  application  of  the  Markov  approximation 


Let  us  consider  corrections  to  the  Markov  approximation 
when  a  finite  radius  of  correlation  in  the  propagating  direct - 
tion  is  taken  into  account.  Let  £(},?)  be  a  sharp  peaked 
function  of  the  argument  J-J'  and  < £  (J,f) id' *  B  (}*}', 

The  parabolic  equation  for  the  random  field  U-(J, f)  is 
written  as: 


2i.*e 


3U(^f) 

*3 


+ 


V?  U  Ci.t)  + 

■  *  *  (A3  *1 ) 


Taking  the  Fourier  transform  of  eq.  (A3.1)  with  respect  to 
and  averaging,  one  obtains 


3  <  u.>  <  £  ( j  (  A  -  ±|?f<  £  (If) 

**•  *  Jei«j 


32 


?)u  (},?)> 

tv 

f 


(A3. 2) 


Applying  the  renormalized  series  expansion  described  in 


Appendix  1,  we  have 

±*±<£  (*,?)«■<$/?)> 

2 


+ - J  J^W(A3  * 3 } 


If  the  correlation  length  in  the  z  direction  is  very  small 
such  that  the  1st  term  in  the  series  expansion  of  eq.  (A2.1) 
dominates,  we  can  approximate  (A3. 3)  as  followings 
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2  ^  ^ 


-  T? j  B  H’-> 

G  f  j.f.  j1/?';  hi<»:*  (A3-4> 

^  »  *■  «> 


One  rewrites  (A3. 2)  as 


,f,  i'.f’Ku  e  * 

Ar  A' 


A* 

•V  A^ 


(A3. 5) 


assigning  the  field  £<3,0  is  statistically  homogeneous , 
so  that  its  two  dimensional  spectral  density  satisfies  the 
relation i 


<£  (3, *,,)>  =  $  (£,+  *f,)F <}-}',  kr, ) 

v  «w  ^  V  1  v 


(A3. 6) 


Applying  (A3. 6)  one  rewrites  (A3. 5)  after  some  manipulations, 
as  follows i 

*°°  £ 

s',  fef(- ?> 


Since  the  function  F(z-z',K^  )  is  a  peaked  function  of  the 
argument  z-z',  one  infers  that  the  main  contribution  to  the 
integral  over  z'  lies  in  the  neighborhood  of  z  *  z’.  We 

A 

shall  expand  the  functions  G  and  <  u >  in  powers  of  2-2%  V121 
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G(3.f,J’,t)  =  Gt  (3,r,$.v  -  <}-i')4$ 


33' 


3*3 


3-3 

+ . 

.A  /»• 


<a(3',^)>=<w.(3,^)>  -(3-3';^|i^|  (A3<8) 

.  (3-3’/  aa<iX(3'  **>)>( 

+  2  5F 


3*3  ^ 


Substituting  (A3. 8)  into  (A3. 7),  one  obtains  after  some  te¬ 
dious  calculations 


=  Ai->c5)  (A3-9) 


In  the  coordinate  representation,  this  equation  takes  the 
form 


—  A1  f  •  A  /■*%  Aj,x(3,)«1  Aa.4^-^ 


■  ■■■  (A3. 10) 


•J  <ul}.t)> 


where 
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Ad(i,r)*  z\l?  <i-}',P)(3-3'f  Jj' 

~  J0 

(A3. 11) 

a 0 

A*.  pd)  =  jJ^(^p;iP|Pdp 

•09 

(A3. 12) 


To  use  the  Markov  approximation,  the  terms  on  the  right  side 
in  eq.  (A3. 10)  have  to  be  small  compared  to  the  corresponding 
terms  on  the  left  in  this  equation.  Thus,  one  can  vrite 
thes  conditions  as  follows i 

to  -fe/A.,.  »  A  a,* 


»)  A 3,  <  <  1 


A  similar  analysis  can  be  performed  for  the  higher  moment 


case. 
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APPENDIX  4 


In  Appendix  2,  we  examined  several  expansions  for  VcG 
and  used  numerical  methods  to  test  their  convergence.  We 
assumed  that  the  imaginary  part  of  G  was  very  small  compared 
to  its  real  part,  i.e.  G  is  approximately  real.  This  is  not 
generally  true.  Hence  all  of  the  results  in  Appendix  2  only 
apply  to  ranges  wherein  G  is  real. 

Behavior  of  G  w.r.t.  Ky(=B+ifc  )t 

We  are  dealing  with  the  simple  case  defined  by  the  equa¬ 
tion  « 

»  A  'V  ^  A 

ijrQ  -  “  1 

(A4.1) 


where  V  is  a  random  number.  Taxing  the  Fourier  transform  of 
(A4.1),  one  obtains  the  algebraic  equation! 


O  +  ifc-V.)?**  1 

(A4 . 2 ) 


1 

B  +  *  «  -  V 


(A4.3) 


Therefore 

A  6-  v 

Real  part  of  *  CB-V^t  *  * 


(A4.4a) 
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a 

Img.  part  of  Ik 


-  6 

(B-v)'i- 


(A4.4b) 


If  the  probability  distribution  function  F(v)  of  V  is  given, 
we  can  evaluate  the  average  Green's  function  exactly i 


and 


<lkr> 


( v )  p(  v*;  d  or* 


< 


*  «o 

y  —  «0 


(A4.5) 


On  the  other  hand,  we  can  expand  Vc^  in  the  following  wayj 

t  jVv*>  +  ^  (A4.6) 

t . 


A 

The  first  approximation  for  4  $  k  ^  is  defined  by  the  follow¬ 
ing  equation i 


(B  +  i*  i 


(A4.7) 


i 

% 


\ 


The  2nd  approximation  is  defined  by: 


vx>  3^-  ( <  v*>  -  x  <  va>  )  3*  J  ^ 

(A4.8) 


etc. 


We  shall  examine  the  approximate  solutions  of  (A4.7)  and 
(A4.8)  and  compare  these  solutions  with  the  exact  solution 
(A4.5)  by  numerical  methods. 

(a)  For  the  case  that  V  is  uniformly  distributed  over  a  range 

from  -a  to  +a,  one  obtains  for  the  2nd  moment  of  V  t 

~:v  __  a.2 

<V  >  »  3 

and  for  the  4th  moment  of  V  t 


Substituting  these  moments  into  eq.  (A4.7)  and  (A4.8)  and 
using  contour  integral  techniques,  we  evaluate  ^  k  as  a 
function  of  B  and  depicted  in  Fig.(A4.1)  to  Fig.(A4.6). 
Explanations  of  Figures: 

Fig  A4.1: 
variance  =  0.333 

dotted  line  *  1st  approx,  of  %  Kr 
solid  line  =  2nd  approx,  of  ^  Jcr 
Fid.  A4.2: 

— -  n  ,  ~x 

The  exact  solution  of  fcr  for  <V>  =  0.333 

Fig.  A4.3  : 

variance  =  2.083 

dotted  line  =  1st  approx,  of  ^  Kr 
solid  line  =  2nd  approx, of  §  kr 
Fig.  A4.4: 

The  exact  solution  of  ^  kr  for  4  V y  =  2.083 


—4  S50«-01  l_ 


-3  aa«-ei  I _ _ _ 

-w  a  -s  aa 


j _ i _ 

a  e  5  aa 


B  _ 

Fig.  A4.5.  g^-  vs.  B;  v  is  uniformly 
over  a  range  (-a,  +a);  <v*>  *  8.333* 
dotted  line;  2nd  approx,  solid  line. 


ia  a 

distributed 
1st  approx. 


v*r«  8  33312 


3  ao«-ai  _ 


Fig.  A4.6.  The  exact  solution  of  g^  for  <v> 
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Fig.  A4.5i 
variance  =  8.333 

dotted  line  =  1st  approx. of  %  kr 
solid  line  ■  2nd  approx,  of  ^  kr 
Fig.  A4.6i 

The  exact  solution  of  $  kr  for  4v>  =  8.333 


Conclusion! 

For  a  r.v.v”  with  a  uniform  distribution  the  approxima¬ 
tion  procedure  displays  quite  good  convergence. 

(b)  For  the  case  that  V  obeys  a  Gaussian  statistics,  one  ob¬ 
tains  for  the  even  moments  of  V» 

^  iw  ,  vn 

i  V  >=  1  3  5  . . .  (2n-l )  <V  / 


Substituting  these  moments  into  (A4.7),  (A4.8)  and  apply¬ 
ing  the  same  technique  as  in  (a),  we  can  evaluate  ^  k  as 
a  function  of  B.  The  results  are  shown  in  Fig.  (A4.7). 
Explanation  of  Fig.  (A4.7)i 
variance  =  1 


dotted  line  =  1st  approx  of 


kr 


crossed  line  (x)  =  2nd  approx  of  $  kr 


kr 


solid  line  =  3rd  approx  of 
thin  line  =  exact  solution  of  ^  kr 
Conclusion! 

We  observe  that  these  approx,  procedures  do  not  show 
convergence.  However,  we  shall  consider  some  other  obser¬ 
vation  before  making  any  definite  conclusions. 

Since  Vc  can  be  evaluated  exactly  from  the  following  inte¬ 
gral! 

A 


(A4.9) 
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We  shall  compare  the  real  and  img.  parts  of  VcG  for  the  1st, 
2nd  and  3rd  approximations,  respectively,  with  the  exact  re¬ 
sult. 

For  the  case  that  V  obeys  a  Gaussian  statistics,  we  display 
a  few  plots  of  these  results i  (  write  g^  as  G  } 

Fiq.A4.8i (exact  result) 
variance  =  1 

curved )  =  real  part  of  VcG 
curve(2)  *  img  part  of  VcG 
curve (3)  =  absolute  value  of  VcG 
Fig.  A4 . 9t 
variance  =  1 

curve  (1)  =  real  part  of  VcG 
curve  (2)  =  img  part  of  VcG 
curve  (3)  ■  absolute  value  of  VcG 
Fig.  A4.10i 
variance  =  1 

curve  (1)  =  real  part  of  VcG 
curve (2)  =  img  part  of  VcG 

curve  (3)  =  absolute  value  of  VcG 
Fig.  A4.ll: 

For  variance  =  0.25 

curve  (1)  -  real  part  of  VcG 

curve  (2)  *  img  part  of  VcG 

curve (3)  *  absolute  value  of  VcG 

Fig.  A4.12i 

For^V1^  •  0.25 

curve  (  1)  *  real  part  of  VcG 

curve  (2)  =  img.  part  of  VcG 

curve  (3)  =  absolute  value  of  VcG 

Fig.  A4 . 13 1 

for  ^»V*>  ■  0.25 

curve  (1)  *  real  part  of  VcG 

curve  (2)  ®  img.  part  of  VcG 


r 
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curve  (3)  =  absolute  value  of  VcG 

The  above  observations  show  that  one  can  not  get  a  good 
approx,  for  VcG  by  taking  a  finite  number  of  terms  in  eq. 
(A4.6)  for  a  Gaussian  statistics.  We  conclude  that  the 
"moment"  expansion  for  VcG  is  not  suitable  for  a  Gaussian 
statistics. 

For  the  case  that  V  is  uniformly  distributed  over  the  range 
from  -0.5  to  +0.5,  we  plot  VcG  both  exactly  and  for  several 
approximations,  viz* 

Fig.  A4.14i  (exact  result) 
variance  =  8.333  x  10“  2 
curve  (1)  *  real  part  of  VcG 
curve  (2)  »  img  part  of  VcG 
curve  (3)  =  absolute  value  of  VcG 
Fiq.A4. 15i  (1st  approx.) 
variance  =  8.333  x  10“  2 
curve  (1)  ~  real  part  of  VcG 
curve  (2)  =  img  part  of  VcG 
curve  (3)  »  absolute  value  of  VcG 
Fig.  A4.16  t  (2nd  approx) 
variance  =  8.333  x  10  “2 
curve  (1)  =  real  part  of  VcG 
curve  (2)  =  img  part  of  VcG 
curve  (3)  =  absolute  value  of  VcG 

We  observe  that  eq.  (A4.6)  provides  a  rapidly  convergent 

A* 

series  expansion  for  VcG  when  V  obeys  a  uniform  statistics. 


-3.ee 


.ee  -i.se  e.e  i.se  see 

®  <%*  vl 

A4.7.  g^.  vs.  B;  v  is  Gaussian;  <v>  =  1. 

..  1st  approx.;  (2)  xxx  2nd  approx. » 

3rd  approx.;  (4)  thin  linei  exact  solution. 


S  -1  V 


Fig.  A4.9.  1st  approx,  of  fig.  A4.8*. 


Fig.  A4.ll.  v  is  Gaussian;  ;  (1)  real  part 
of  Vcg;  (2)  Img.  part  of  V_.g;  (3)  absolute  value 
of  Vcg.  <  v\>  =  0.25. 


Fig.  A4.12.  1st  approx,  of  Fig.  A4.ll. 


3.i 


•Fig.  A4.13.  2nd  approx,  of  Fig.  A4.ll. 


3.1M-66 


200  - 


APPENDIX  5 


Limits  of  Application  of  The  Parabolic  Equation 

The  parabolic  equation  differs  from  the  complete  scalar  wave 
equation  by  a  term  5*u/®3*  •  one  treats  u/  as  a 
small  perturbation  term  and  writes 


2  i-fe. 


A 

*3 


3  A  .  Z  ~  A 

+  7f  = 


3  3" 


(A5.1) 


then  one  expands 


u<i,n=  u, a.f)  +  u»(5'5’->+ .  <A5.2) 


and  supposes  that  u^  is  of  the  same  order  of  smallness  as 


Substituting  (A5.2)  into  (A5.1)  and  equating  groups  of 
terms  of  the  same  order,  one  obtains 

A 


2ifc.Lii.  +  +  fc.1  £  w- i  =  o 


J  V  ^ 


3  U, 


2-U„|4i  +  +  *•  e  u‘  =  '  aj 

*  5 


(A5.3) 
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with  the  boundary  conditions 


u,  («,  n  *  u„(f)  and  U,(o,n  =  O 

*  V  V  *  ^ 


one  finds  on  taking  the  ensemble  average  of  (A5.3) 


+  ^  +  £  u,>  ~  o 


(A5,4a) 


3  V  U,  > 

*y 


(A5.4b) 


Using  the  renormalization  procedure  described  in  Appendix  1> 
one  writes  equation  (A5.4a)  as 


(l  ,-vc)<Cl>  =  o 

on  using  of  the  approximation 


'  .  y.  , 


vcu  S'  J 

u  a'.?1)  di'jy ' 

<%»» 

and  the  assumption  that 

<£0 ,f>e  (}'?'»  =  £a-3')A(t-i'> 


(A5.5) 


(A5.6) 
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One  obtains,  after  some  manipulation, 

vc  u  =  -  u  (3,?; 

4  ~  2 

s  -  ihl  A  ( 0 -> 

S' 

Therefore,  we  can  rewrite  eq.  (A5.4a)  as  follows: 

+  vr*L +  ^  w>  =  ° 

Similarly,  eq.  (A5.4b)  can  be  rewritten  as 

+  Vf*u,  +  =  -fjr 

Let  us  consider  the  case  of  a  plane  wave  incident  upon 
dom  medium  such  that 

U0(?)=  Can.stvsMi't  =»  Ut» 

One  obtains,  from  equation  (AS. 4a),  that 

_  <LM(P  ? 

U,  (  l)  a  Ue  €  S 

-  -LotS 

=  u.  e  2  * 

<=<  =  #.*  *  (•> 


(A5.7) 


(A5.8 ) 


(A5.9) 


a  ran- 


(A5.10) 


where 
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equation  (A5.4b)  then  takes  form 


The  solution  of  eq.  (A5.ll)  corresponding  to  the  boundary  con¬ 
dition  u2(0)  =  0  can  be  calculated  as  follows: 
set 


U 


Z 


(A5.12) 


One  obtains  from  eq.  (A5.5) 


U2  (o)  +  (D) 


(A5.13) 


Therefore 


=  u.e 


128 

u'  TIs  <*■ 


(A5.15) 


In  this  case 


U(5>  =  u,  (3>  ■+  + . 


(A5.16) 
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The  field  u, should  be  considered  only  in  the  region 
where  it  is  marking  different  from  zero,  i.e.  for  o<  1»  thus, 
we  obtain  the  condition  for  a  small  correction  term, 

-£.4-  1 

;  ft*  c* 


oL 


Xd  <<  1 


(A5.17) 


If  condition  (A5.17)  is  valid  then  the  correction  term  in  eq. 
(A5.16)  will  be  small  compared  to  its  first  term.  For  the 
higher  moment  case,  the  calculation  is  quite  complicate.  For 
example,  in  the  2nd  moment  case,  one  considers  the  corrections 
to  m2(  Z,  )  =<£  (  2>  f,)  $*(  2 If  u  *  u4  +  U.+  .  .  .  , 

one  has 


<  w,  c$>  *?>>  + 

^  ?i)  u.,  ^  )  y  +  — 

«%*• 

=  m2  +  m2  +  m'2>  + - 


(A5.18) 


Applying  the  Markov  approx,  and  some  very  tedious  calculations, 
one  can  obtain  equation  for  M2  »  M2  ,  M2  . .  as  followsj 


(A5.19) 


0 


{ 2  i  ^  +  ( Vfil-  vf;j)  t  i*t  [A  («;-  A  <  f.  ■ -«] j 


If  the  parabolic  approximation  is  used,  the  correction  terms 
and  M*^*  have  to  be  small  compared  to  the  1st  term  . 
One  obtains  from  eq,  (A5.20)  and  (A5.21),  let  *  M2*' 

the  following  equation  for  M2  » 


Let  us  now  consider  the  case  of  a  plane  incident  wave,  one 
obtains 


-  &{ac*)-A  (r.-faJZ 

i«.*i  e  +  L  ~~ 

(A5.23) 


In  this  case 

^  ^  V 


(A5.24) 


Moreover  • 

(  V-  yn)M-'  =  « 

(  <-  vf;;  mT  =o 
(  vf;  +  vfl2;  m;°  =  2  ^  M? 


Equation  (A5.22)  thus  takes  the  form 


H(S)J 


4 


Vr*M 


C) 

2 


(A5.25) 


(A5.26) 


from  (A5.23),  one  finds 
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M'”'  =  m"’  AJO  I  53  Cn'(f>]2 
2  8  ‘  4? 


_  **3*yf»H  <?;  ? 

*  ) 


(A5.27) 


Therefore 


A(o7 

s ■ 


2_ 


4o$:v^HCf)  [ 
2 


(A5.28) 


this  leads  to  the  conditionj 


if 


then 


At*)  j  VtfH '(?)! 

8  1  4? 


M 


<< 
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Appendix  6 


The  equations  to  be  solved  arei 


(A6.1) 


-  — 1 _ L_  fv^  _  "T  £  M. 


a  r,  s  n 


(A6.2) 


One  writes 


M«  ( Jt*  JJ-  M,  <  f-  -£r  M;  -  rf  M„] 


dr,  3 r, 


(A6.3) 


+  [f* (-  ITT^  M*  -  rf 

-  ,  -%*  (A6.4) 

. 
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The  difference  equations,  to  the  2nd  order  accurracy  in 
are  i 


M«  +  -rf  M„] 

+  [-  rfMxj- 

<rf  [-  Mi  -  *rf  MH]]£l£ 


(A6.5) 


and 


=  Mit  { dn"d r* 

i -iU- 

(rf*$)  H 


(A6.6) 


*'J) 


where 
i 


dr.a  forward  difference  operator  w.r.t.  r. 

*  «• 

d"  =  baclcward  difference  operator  w.r.t.  r. 
*  * 
<*r.  *  center  difference  operator  w.r.t.  r». 

A  * 


The  boundary  conditions  for  M4  are  described  in  equations 
(4.46)  and  (4.47). 
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